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ABSTRACT 


Investigation of problems related to the design and control 
of a mobile planetary vehicle to implement a systematic plan for 
the exploration of Mars has been undertaken. Problem areas 
receiving attention include: vehicle configuration, control, 

dynamics, systems and propulsion; systems analysis; terrain model- 
ing and path selection; and chemical analysis of specimens. The 
following specific tasks have been under study: vehicle model 

design, mathematical model of vehicle dynamics, experimental vehicle 
dynamics, obstacle negotiation, electromechanical controls, remote 
control, collapsibility and deployment, construction of a wheel 
tester, wheel analysis, payload design, system design optimization, 
effect of design assumptions, accessory optimal design, on-board 
computer subsystem, laser range measurement, discrete obstacle de- 
tection, obstacle detection systems, terrain modeling, path selec- 
tion system simulation and evaluation, gas chromatograph/mass 
spectrometer system concepts, chromatograph model evaluation and 
improvement. 

These tasks are defined in detail and the progress which has 
been achieved during the period July 1, 1972 to June 30, 1973 is 
summarized. 
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Analysis and Design of a Capsule landing System 
and Surface Vehicle Control System for Mars Exploration 


I. Introduction 

Current national goals in space exploration include a detailed study 
and examination of the planet Mars. The effectiveness of Mars exploration 
missions would be enhanced according to the extent to which the investigative 
devices which are landed are mobile, to the range of their mobility, and to 
the ability to control their motion. In order to achieve basic mission objec- 
tives, and beyond that, to maximize the return on the commitment of resources 
to the mission, formidable technical problems must be resolved. The major 
factor contributing to these problems is the round trip communications time 
delay between martian and earth control stations which varies from a minimum 
of about 9 minutes to a maximum of 25 minutes. This time delay imposes strin- 
gent requirements on the vehicle, on its control and communication systems and 
on those systems included on board to make the scientific measurements, in 
terms of theii ability to function autonomously. These systems must be able 
to operate with a high degree of reliability and must be capable of calling 
for earth control under appropriate circumstances. 

A number of Important problems originating with these factors and 
relating directly the basic mission objectives of an unmanned exploration of 
Mars have been and are currently being investigated by a faculty- student pro- 
ject team at Rensselaer Polytechnic Institute with the support of NASA 
NGL-33-018-091. 

This progress report describes the tasks which have been undertaken 
and documents the progress which has been achieved in the interval July 1, 

1972 to June 30, 1973. 

II. Definition of Tasks 

The delay time in round trip communication between Mars and Earth 
gives rise to unique problems relevant to martian and/or other planetary 
explorations. All phases of the mission from landing, the capsule in the 
neighborhood of a desired position to the systematic traversing of the sur- 
face and the attendant detection, measurement, and analytical operations 
must be consunmated with a minimum of control and instruction by earth based 
units. The delay time requires that on board systems capable of making 
rational decisions be developed and that suitable precautions be taken against 
potential catastrophic failures. Four major task areas, which are in turn 
divided into appropriate sub-tasks, have been defined and are listed below. 

A. Vehicle Configuration. Control. Dynamics. Systems and Propulsion . 

The objective of this task is the design of a roving vehicle for 
the exploration of Mars. This design includes the aspects of 
vehicle configuration, collapsibllity for launch configuration, 
deployment, dynamics, motion and attitude control, capability 
definition, and performance evaluation. The efforts of this task 
are made within the text of the mission definition and delivery 
system constraints, *he ultimate goal 1 i the development of a 
vehicle capable of operating within the constraints of the Ma-s 
mission, but flexible in design, to Insure reliability with respect 
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to the unknowns of Mars, and to accommodate the alterations 
in the mission as a result of information gathered during the 
exploration. 

B. General System** A naly sis. The objective of this task is to 
develop a framework within which decisions in design involving 
conflicting requirements can be made rationally and in the con- 
text of the whole system and mission. Relationships between 
alternative mission profiles and specifications and weight, 
energy and space allocation and management will be sought. 

C. Surface Navigation and Path Control . Once the capsule is landed 
and the roving vehicle is in an operational state, it is neces- 
sary that the vehicle can be directed to proceed under remote 
control from the landing site to specified positions on the 
martian surface.. This task is concerned with the problems of 
terrain sensing, obstacle detection, terrain modeling, path se- 
lection and navigation between the initial and terminal sites 
when major terrain features precluding direct paths are to be 
anticipated. On board decision making capability must be designed 
to minimize earth control responsibility except in the most adverse 
circumstances. 

D. Chemical Analysis of Specimens . A major objective of martian 
surface exploration will be to obtain chemical, biochemical or 
biological information. Many experiments proposed for the 
mission require a general duty, gas chromotograph-mass spectro- 
meter for chemical analysis. The objective of this task is to 
generate fundamental data and concepts required to optimize this 
chemical analysis system. 

III. Summary of Results 

Task A. Vehicle Configuration. Control, Dynamics, Systevas and Propulsion 

This broad task has been divided into the following subtasks: vehicle 
model design, construction and evaluation; mathematical modeling of vehicle dyna- 
mics; experimental dynamics analysis; obstacle negotiation; electromechanical 
control systems and remote control; collapsibillty and deployment; wheel testing 
apparatus; and wheel analysis. A brief summary of the progress achieved in the 
prior interval in each of these areas follows below. 

A. I. Vehicle Design 

The fundamental features of a "third-generation design" vehicle 
have been developed and are shown to meet collapsibillty, capsule 
limitation, and launch configuration requirements. A deployment scheme 
which is reliable and requires no special terrain for landing has been 
developed. The present design allows for an emergency maneuver in which 
the vehicle can, in fact, reverse itself and provide the means for the 
vehicle to extricate itself from an impossible situation. A conservative 
estimate of the payload fraction of this vehicle is shown to be equiva- 
lent to that estimate for a six wheel vehicle. A 0.2 scale vehicle has 
been built to demonstrate these capabilities. 

A.l.a. Design. Construction end Evaluation of 0.4 Scale Model 


Construction of a 0.4 scale model of the RPI-MRV has been 
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essentially completed. The proposed deployment concept has been 
Implemented and shown to be effective. This model has been used 
for preliminary dynamic studies, the results of which are expected 
to modify the design detail. ft is intended that this vehicle will 
be fully automated and be tested dynamically on appropriate terrains 
by remote control during the coming year. 

A.l.b. Vehicle Dynamics 

Equations describing the dynamic behavior of the RPI-MRV using 
the bicycle configuration have been developed and are in the process 
of being solved by digital computation. In addition, an analytical 
model describing a three-dimensional three degree of freedom model is 
under development. The 0.4 model has been tested on an electromagnetic 
shaker and the experimental results which have been obtained are being 
used to modify the detailed design of the vehicle. These mathematical 
studies combined with the experimental studies are intended to provide 
a sound basis for modifying the 0.4 model and at the same time permit- 
ting extension of the design to a prototype. 

A.l.c. Obstacle Negotiation 

A digital computer program has been developed to model a two 
dimensional vehicle operating in a vertical plane as it traverses 
steps, slopes and crevasses. The program provides knowledge on energy 
requirements to overcome frictional and static torques. In addition, 
the program includes a stability factor which indicates under what 
obstacle circumstances the vehicle's stability may become marginal. 
These calculations and others to be made in the future will provide a 
basis for sizing the motors for the vehicle. 

A.l.d. Radio Control 


The objective of this task is to design, construct and implement 
a remote control capability for the 0.4 scale model. The manual input 
device, command encoder, radio link and command decoder have been de- 
signed, built and tested with the exception of the clocks required at 
the transmitter and receiver to synchronize the data reception with 
the data transmission. Progress has been made on the control system 
in the areas of speed control, command storage, and conversion to 
analog levels. To complete the systems, sensors must be placed in the 
vehicle and their output incorporated into the control systems. The 
radio link and motor controls systems have been Interfaced and used 
to control the vehicle motors at three output speed levels in both the 
forward and reverse directions. As a first step, it is planned that the 
complete radio system be installed to permit open loop manual control 
of the vehicle. This will be expanded first to close loop control and 
ultimately to computer operation. 

A. 2. Collapsibility and Deployment 

A solution for the problem of stowage of the RPI-MRV in the 
Viking capsule and its subsequent deployment on Mars has been de- 
veloped. The vehicle is folded in such a manner as to locate the 
center of the payload fraction on the centerline of the capsule, and 
to locate the wheels in appropriate positions. The deployment concept 
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is such as to permit deployment on relatively difficult terrains. 
Furthermore, the particular configuration which has been developed 
allows the vehicle to reverse itself, if necessary, to avoid an im- 
possible obstacle situation. 

A. 3. a. Wheel 'fester and Grouser Design 

A wheel tester system has been designed, built and used to test the 
RPI toroidal wheel and alternative grouser designs. Slip and power tests 
have been performed and evaluated for five alternative grouser designs. 

A . 3 . b . Wheel Analysis 

Existing mathematical models for the design of the toroidal wheel 
have been appraised and evaluated. A new iterative mathematical analy- 
sis has been developed and has been applied to experimental data. Com- 
parison with experimental data and predicted results for concentrated 
loads have proved to be favorable. Future work will be directed to the 
problem of distributed loads. 

Task B. General Systems Analysis 

This major task is subdivided into four subtasks: system modeling and 

design optimization; sensitivity of optimal designs; time simulation of vehicle 

performance; and specification of an integrated MRV controller. 

B. l. System Modeling and Design Optimization 

During the past year, the objectives set for this subtask have 
been achieved. Optimal designs have been identified and the implica- 
tions of these results have been investigated for three major alterna- 
tive models, i.e. four-wheeled vehicle with direct earth link, six- 
wheeled vehicle with direct earth link and a four-wheeled vehicle with 
an orbiter-earth link. These optimal designs specify values of the 
parameters of the system, e.g., weifit am* power allocations to the 
vehicle subsystems and subsystem performance levels (communications 
data rate, vehicle velocity, science time, etc.). It is believed that 
the overall approach to the system analysis problem which has been 
developed represents a contribution to the design methods for large 
aerospace systems. 

B.2. Sensitivity of Optimal Designs 

The method for the linear approximation of effects of forced 
perturbation of one of the parameters of a design from its optimal 
value has been derived. The sensitivity of the optimal rover designs 
which have been obtained are being investigated by this method. 

B.3. Time Simulation of Vehicle Performance 

A simulation program that allows the investigation of the per- 
formance of designs determined under B.l. in real-time has been 
written. It is planned that these optimal designs be studied in the 
real-time sense to determine their feasibility. The objective is to 
analyse the performance of the designs identified under B.l. with 
differing mission profiles end stochastic inputs, (terrain slopes, 
emergencies, etc.). 
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B.4. Specification of an Integrated HKV Controller 

The objective of this subtask is to investigate the performance 
requirements ami candidate designs for the MRV on-hoard computer. A 
queuing theory structure lias been developed for determining the type 
of computer architecture appropriate for a given set of requirements. 

Task C. Terrain Sensing, Terrain Modeling, Obstacle Detection and Path Selection 

This task deals with problems of terrain sensing, terrain modeling, 
obstacle detection and path selection for an autonomous roving vehicle. Specific 
tasks receiving active effort include: the design of a laser rangefinder for 
terrain sensing; design of a laser scan system; stochastic estimation of gradient; 
stereo and non-stcreo obstacle detection systems; and path selection systems 
evaluation. 


C.l.a. Laser Rangefinder Design: Pulse Amplification 

Various systems by which a laser may be used for rangefinding 
purposes in the range of 3-30 meters have been analyzed and evaluated. 
Optical focusing, timed pulse and phase comparison techniques were 
compared. It has been concluded that a pulse amplification technique 
or a hybrid system combining the timed pulse and phase comparison con- 
cepts are most promising. Of these two systems, pulse amplification 
appear to be more suited for rapid terrain scanning. A preliminary de- 
sign of a rangefinder based on pulse amplification has been completed 
and some of the required electronic components have been constructed 
and tested. 

C.l.b. Laser Scan Systems and Laser Rangefinder Design: Phase Comparison 

Alternative laser scan systems have been considered and an 
alternative utilizing digital deflection optics or a matrix of laser 
diodes and a vibrating detector mirror which can meet the rapid scan 
requirement has been developed. This scan system can support a pulse 
amplification rangefinder system su h as that described under C.l.a. 
above, or a system based on phase comparison as proposed under this 
task. 


A design of a laser rangefinder based on phase comparison has 
been developed and construction and testing of some of the electronic 
components has been completed. Unfortunately critical components could 
not be constructed during this past period and a decisive appraisal of 
the feasibility of the phase comparison method will have to be defer- 
red to the future, 

C.2. Terrain Modeling end Obstacle Detection 

C.2.a. Error Analveie end Stochastic Pit for Gradients 

An analysis of the effect of measurement errors, i.e. , range 
end elevation and azimuthal angles, on the gradient of terrain in 
the 3-30 meter range has been undertaken. It ie found that measure- 
ment errors originating with vehicle motion, (l,e. roll, pitch end 
heave) result in very large errors in estimating the in-path end 
cross-path slopes and the gradient. It ia concluded that e rapid 


scan system should be used to eliminate those errors. Assuming a 
rapid scan system obviating vehicle motion effects, the rein 1 
ships between the standard deviation of eradiemt as a funci^-u of 
range and of the standard deviation range and elevation am azimuths, 
error for a horizontal flat plain have been detei mined. The standard 
deviation of gradient varies from less than 1° to about 8° depending 
on the measurement errors of up to a standard deviation of 10 cm "n 
range and for a data spacing corresponding to elevation and azimuthal 
angle increments of 1 minute. Data spacing as determined by elevation 
and azimuthal angle increments has a major effect on the gradient 
error. As the spacing is increased, the gradient error is reduced. 
Unfortunately large data spacing also implies increased loss of infor- 
mation. A trade-off between data spacing and gradient error will be 
necessary. The effect of an inclined terrain has also been deter- 
mined. 

As this effort more completely relates the effect of measurement 
errors on the reliability of terrain model interpretation, specifica- 
tion for a terrain sensing system can be defined more rationally. 

C.2.b. Gradient Estimates Based on Stereo Measurements 


The method of calculating a stochastic estimate of the gradient 
has been applied to two stereo-measurement systems and compared to 
the non-stereo range-angle system. Instrumentation accuracy require- 
ments to meet a 2 degree gradient standard deviation have been deter- 
mined for each alternative. The non-stereo system is shown to be 
superior to both the stereo-range and stereo-angles alternatives and 
indeed is the only system which can be conceivably implemented with 
only modest improvements in the state of the art of laser range 
measurement. 

C.3. Path- Select ion System Evaluation 

The principal efforts have been devoted to improving the 
realism of the digital computer simwatlo*. program developed in the 
previous year and in using the program for system evaluation purposes* 
Specifically, simulation of an emergency backup mode, non-zero vehicle 
dimensions, and indices for safety and obstacle detection errors have 
been incorporated into the program* In addition, the program has been 
used to compare two different mounting arrangements for the range- 
measuring sensor used by the path- select ion system* The two configura- 
tions compared have the plane of the sensor beam fixed relative to the 
vehicle and fixed relative to the local planet vertical* 

Task D. Chemical Analysis of Specimens 

This task is concerned with developing the fundamental concepts which 
will be required to optimize a gas chromatograph-mass spectrometer analysis 
system* Topics receiving active effort include mss spectrometer characteristics, 
carrier gas gene ratio;: and removal, multicomponent chroMtography , and chroma- 
tograph model improvement* 


D.l 


Gas Chromatograph/Mass Spectrometer System Concepts 


7 . 


D. 1 , a . Mass Spectrometer System Char a cter 1st ier. 

The mass spectrometer studies are intended to characterize die 
major properties of die instrument, the resulting relations bein;; 
suitable for investigation of system power, weight, and volume re- 
quirements. Relations involving mass resolution for single and 
double focussing instruments have Leen developed. Preliminary sLudies 
showed the theoretical analysis of instrument sensitivity to be beyond 
the scope of this project, and empirical relations are being sought. 
Investigation of the effect of major design variables upon power re- 
quirements ras begun. Preliminary analysis of the characteristics of 
the ionization chamber filament was completed. 

D. l.b. Carrier 'las Generation and Removal 


Preliminary analysis of high-pressure gas storage as a source 
of carrier gas has been completed. Use of filament-reinforced con- 
struction will provide storage capacity for one year's operation of 
the chromatograph within a volume of about one or two liters and a 
weight of about 2 kgm. The concept appears competitive with a pre- 
viously proposed electrolytic gas generator. Engineering analysis of 
carrier gas removal by diffusion through a palladium alloy tubing be- 
gan. A system containing only carrier gas and no sample was successfully 
modeled. Future work with the model will incorporate the effect of 
the nondiffusing sample. 

D.2, Chromatographic Systems Analysis 

D.2.a. Multicomponent Chromatography 

Data for the system n-pentane/n-heptane on porous Chrotnosorb-102 
adsorbent were obtained at 150, 175, and 200° C for mixtures contain- 
ing zero to 100% n-pentane by weight. Prior results showing limita- 
tions on superposition of pure component data to predict multicomponent 
chromatograms were verified. A nonporous adsorbent failed to separate 
the system because of large input sample dispersions. Modifications to 
the system to improve sample injection are being considered. An auto- 
mated data processing scheme involving magnetic tape recording of the 
detector signals and processing by a minicomputer was rejected because 
of resolution limitations of the available A/D converters. Preliminary 
studies on porosity and pore size distribution of the adsorbents were 
begun. 

D.2.b. Chromatograph Model Improvement 

Intraparticle diffusion and rates of adsorption were Incorporated 
into the mathematical model of the chromatograph to improve representa- 
tion of the experimental data. Moment analysis of the resulting system 
of equations showed the model to be capable of representing store ade- 
quately the spreading of experimental data. It appeared an analytical 
solution of the model was not realizable, and numerical solution methods 
are being investigated. Criteria for substituting a finite boundary 
condition required for numerical solution in place of an infinite bound- 
ary condition used in prior analytical work were developed. Computational 
schemes based on finite differences ware shown to be inefficient, so the 
method of weighted residuals with orthogonal collocation is being investi- 
gated. 
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IV. Detailed Summaries of Progress 

Task A, Veln'clc Conf iguruLion, Control, Dynamic.-., System and Propulsion 

This broad task is concerned with all aspects of the design of a 
roving vehicle for Mars exploration. Included within this major effort are 
the following subtasks: vehicle model design, construction and evaluations; 
mathematical model of the vel '-.le dynamics; experimental dynamic analysis; 
obstacle negotiation; electromechanical control systems and remote control 
systems; collapsibility and deployment; wheel testing apparatus; wheeL analysis; 
and payload design. 

Progress made in these areas during the past interval is summarized 
in the sections which follow. 

A . 1 . Vehicle Design 

Student Group Leader: J. Almstead 

Faculty Advisor: Prof. G. N. Sand' 

The objective of this task is to develop the concepts of the RPI- 
MRV ar.d evaluate these concepts re 1 -_ive to the rover mission require- 
ments. Efforts during the past year have been directed toward making 
the fundamental features of this "third-generation design" meet the 
collapsibility, capsule space limitation and launch configuration re- 
quirements, developing a deployment scheme which is reliable and requires 
no special landing terrain, and refining the emergency maneuver to better 
protect the instrumentation when executed. While satisfying these three 
areas of requirements, the design had to remain as simple as possible to 
keep structure weight minimized and reliability maximized. 

Figures 1 through 8 show the vehicle concept which was developed. 

It has the capability of collapsing to fit in the Viking capsule. When 
in the capsu'.e, the payload is centered and attached to the mounting 
ring in the capsule bottom. The frame is locked in place and braced 
to the capsv'e walls. The wheels are compressed slightly to act as mounts 
in the closed capsule, and provide a positive force for the separation of 
the capsule. 

Deployment on Mars requires only three vehicle motions which can be 
done as long as the vehicle has landed safely and the payload is stable 
on the ground. Because the wheels are raised above the ground, deployed, 
and then placed on the ground, the limitations on the terrain at the 
landing site are minimized. The three vehicle motions Involved in the 
deployment are rear strut deployment, frame unfolding, and front axle 
swing-over. The details of the deployment were reported in Ref. 1. 

The emergency maneuver, also described in Ref. 1, permits the 
RPI-MRV to escape from situations which would otherwise terminate the 
mission. The vehicle is now able to execute the maneuver without having 
to use the payload as a skid as did the second generation version. With 
symmetry, the vehicle is able to switch front for rear and rear for front, 
and drive away from the danger. 

A first approximation of the payload fraction was made. It defined 
payload fraction to bet (Total permitted weight-vehicle) /(Total per- 
mitted weight), or effectively: (payload weight)/ (total permitted weight). 
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Fig. 1. The RPI-MRV in collapsed launch/land configuration. 
This 0.184 scale model fits into a similar scale 
model of the Viking capsule. 
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Figure 3. With payload package still on the ground, and the 
front structure and wheels deployed, the rear wheels 
begin to unfold. 
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Figure 6. With the payload still on the ground and the rear 

wheels rotated, the front frame and wheels are swung 
to the opposite side of the payload, clearing all 
antennas and vertical extensions of the payload. 
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Figure 7. The front wheels touch ground 'n the rear of the 

I ' payload. Note the payload configuration relative to 

{ the new front and rear wheel locations. 
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Figure 8. With all four wheels again on the surface, the 

payload is raised to desired height, and the RPI-MRV 
drives away from the obstacle. 
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Included for the calculation of the vehicle weight are estimates of 
the frame weight, motor weight and frame link weight. The results in- 
dicated that the greatest contribution to the vehicle weight was due 
ti the motors. Those weights were conservative estimates, which in- 
cluded the motor, motor mounting and gearing weights. The motor weight 
estimate was 65% of the total vehicle weight estimate. The calculated 
payload fraction was 0.585. The indication is that even with the very 
conservative motor weight estimates, the payload fraction of the RPI- 
MRV is equivalent to the payload fraction estimation made for the six 
wheel vehicle, Ref. 2. 

A 0.2 scale vehicle and capsule model has been built to demonstrate 
the vehicle capabilities. The basic design concepts are well established 
for the RPI-MRV. What remains is to evaluate the vehicle design rela- 
tive tc the mission requirements. 

In the finalizing of the vehicle design, the areas of concern are: 
evaluation and testing of existing mechanisms, design of locking and 
deployment actuation mechanisms, physical damper designs and investiga- 
tion of operating parameters, such as power and torque requirements, 
structural components, payload fraction and science package operation. 


A, l.a. Design, Construction, and Evaluation of 0,4 Scale Model - W. Embleton 

Faculty Advisor: Prof. G. N. Sandor 

Construction of a 0.4 scale model of the RPI-MRV has reached a 
stage which will accommodate the first tests and evaluations of the 
larger scale vehicle. This modelling has indicated the need for addi- 
tional design work in the areas of lateral stiffness of the wheel, front 
suspension and motor mounts. It has also provided information on the 
constraints associated with the hinges, joints and locking mechanisms 
for collapsibility and deployment. While bringing to light problems in 
the areas mentioned, the scaling up to th^ 0.4 scale has reinforced the 
RPI-MRV concept as a plausible candidate for the rover mission. 

The requirements for the structure of the 0.4 model were deter- 
mined by area-moment analysis. Materials were selected to minimize 
the frame weight. Once constructed, it was found that the section used 
for the frame was inadequate because of its response to shaker inputs, 
contributing to the oscillations of the vehicle. Also, the frame pro- 
vided insufficient support for the frame pivot. It was changed, the new 
frame being more rigid with the frame pivots removed to accommodate the 
shaker testing. This frame is temporary. When complete, the 0.4 scale 
model will have the required pivots and locking mechanisms for implementa- 
tion of the collapsibility and deployment scheme as well as the emergency 
maneuver. 

The problem of the lateral softness of the wheel was noted during 
the shaker tests. The vertical input to the front was resulting in a 
significant horizontal response of the rear of the vehicle. An investi- 
gation is presently underway which will evaluate the lateral softness 
problem and make recommendations as to the most effective approach to 
increasing lateral stiffness of the wheel while retaining its radial 
compliance and minimal weight. 
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The front suspension originally on the 0.4 scale model proved to 
be too soft when the vehicle was tested on the shaker unit. It was 
replaced with a rigid axle. The shaker data taken in both conditions 
will be compared to indicate the effects of a soft axle versus a stiff 
axle on vehicle response. 

The motor mounts and motor locations will depend on the kind of 
motors used to control the vehicle functions. Presently, an effort is 
being made to locate appropriate DC motors for the 0.4 model. 

Future work will be directed to designing, constructing and evalua- 
ting frame pivots and locking and deployment actuators and to providing 
the propulsive power and control systems required to evaluate the per- 
formance of the vehicle. Efforts to optimize the design will continue. 

A. l.b. Vehicle Dynamics - Mark C. Rodamaker, Alan Ryder 

Faculty Advisor: Prof. G, N. Sandor 

The objective of this task is to provide a mathematical model and 
computer simulation of the vehicle to: 1) aid in the basic design of the 
full scale vehicle, and 2) optimize the full scale design by predicting 
the dynamic response of the vehicle to various inputs. Two approaches to 
the development of the analytical model are currently being studied. Both 
approaches are described below. 

The first analytical approach assume the vehicle to be of a bicycle 
configuration. Equations for this six-degree-of-freedom planar model 
have been derived using a complementary virtual work method. This allows 
the inputs to be displacements rather than forces as required in classi- 
cal Newtonian methods. An analog simulation of the model was attempted. 
However, technical problems resulting from the mating of more than one 
analog computer forced this method to be abandoned. A digital computer 
package "DYNAMO" (DYNAmic Models) is currently being used in an attempt 
to solve the derived equations. 

The second, more recent analytical model considers a three-dimensional 
three-degree-of-freedom model (see Fig, 9). Equations have been developed 
for a moving reference frame as follows: 


y “ vertical translation 

0 X ■ rotation about longitudinal axis (roll) 

0 f ■ rotation about lateral axis (pitch) 

With these » the equations of motion take the form: 


V - «y(2vw + e x a 2 K rW " V bz r ~ •V 
♦ V 4 A> -.rtVrtV + V-‘i*r *2W> 

1 r« 6 . + * y(bZ r ~»Z,) + $2<-b*Zj-.*Zj4v r h) 
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where: 


Z = rear suspension mechanical Impedance (in the vertical 

direction) 

Z^ *» front suspension mechanical impedance (in the vertical 

direction) 

Kj «* mechanical impedance of side suspension (in the vertical 

direction) (front-rear wheel combination) 

m r ,w r ■ mass, weight of sprung vehicle 

h “ height of center of mass of sprung vehicle above 

horizontal plane of unsprung center of mass 

£j,£, 2 »a»h ■ dimensions locating center of sprung mass, Figure 9. 


These simultaneous equations remain to be linearized and solved. 

Actual physical testing of the 0.4 model has been conducted. A 
capacitance-displacement measurement system has been used to measure 
responses. An electro-magnetic shaker unit has provided oscillatory 
inputs to the model. Three dimensional front end frequency responses 
have been measured, in addition to step response data for the front end. 

Rear frequency responses have been practically determined. Vehicle 
characteristics (mass, spring rates, damping coefficients, nrmients of 
inertia) have been measured. 

This will serve as the background for analyzing the measured dynamic 
responses in order to: 

1) evaluate and refine the analytic models 

2) determine the source of undesirable response characteristics. 

In the future, it is proposed that the three-dimensional mathematical 
dynamic modelling of the 0.4 scale model be extended to include additional 
degrees of freedom. Computations wil. be c-jmpared with experimental re- 
sults to verify the mathematical analysis. The analysis will then be used 
to determine the spring and damping coefficients to be put on the experi- 
mental 0.4 scale model. 

In the development of the bicycle model, the next stage Involves 
modeling of terrain features. This is necessary to calculate the dynamic 
pitch and heave response of the RPI-MRV to "typical" Martian surface inputs. 
This terrain model will later be applied to the extension of the bicycle 
model: a three-dimensional multi-degree-of- freedom dynamic model. Loading, 

calculated from the three-dimensional analysis, will be used to else the 
structural members of the RPI-MRV. Analysis of the structurally sired 
vehicle will permit reliability and payload fraction analysis. 

A.l.c. Obstacle negotiation - Larry Mabius 

Faculty Advisor: Prof. 0. N. Sandor 

The objectives of this task were to extend the obstacle negotiation 
study to include the variability of the vehicle center of gravity, a 
stability factor for the vehicle condition, a generalisation of the 
obstacles to a series of straight lines in a plane, and background provision 
for the eventual extension of the Obstacle negotiation modeling to three 
dimensions. 
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A digital computer program was developed to model a two dimensional 
vehicle operating in a vertical plane as it traversed steps, slopes and 
crevasses. The geometric technique used in earlier studies was abandoned 
in favor of the more general force and moment equation, which simplified 
the analysis and accommodates the inclusion of an adjustable payload 
height as an additional model parameter. The study indicated that by 
lowering the eg. of the vehicle, (at the same time lengthening the vehicle 
and placing a greater percentage of weight on the front wheels') the rear 
wheel drive torque required to drive the front wheels up a simple obstacle 
increased. This was expected, because of the increase in loading of the 
front wheels while lowering the center of gravity. The rear wheel drive 
torque required for simple step, crevasse and slope obstacles decreased 
as the eg. was lowered. The model confirmed what had been expected. This 
means that, if needed, an advantage may be gained by reducing the rear 
wheel torque required to traverse an obstacle by eg. height adjustment. 

The torque requirements calculated are based on static equilibrium, not 
including friction or soil-wheel interaction. 

The program was extended to include a stability factor. The intent 
was to have this factor serve as an indicator of how close the center of 
support (vertical projection of eg.) came for some attitude and configura- 
tion to the wheels (the vehicle support points). This will then give, 
some indication of the magnitude of disturbance the vehicle could tolerate 
before becoming unstable. (Instability is here defined as the vertical 
projection of the eg. being outside the quadrangle of wheel contacts on the 
ground). 

A friction term was also included in the model by indicating a 
torque necessary to overcome the rolling friction of the vehicle. As data 
becomes available, this term will be altered to Increase its realism for 
the vehicle. Presently it is included as a constant, but could be made a 
variable dependent upon vehicle life, speed, soil conditions, and wind 
velocities and direction. 

Energy requirement calculations were made based on the frictional 
and static t rques. This should be extended to energy calculations which 
include the effects of factors of terrain, soil conditions (including 
slip), relative wind velocities, and vehicle speed. 

Additional work in this area will be aimed toward developing an 
extension of the planar two-dimensional vehicle and obstacle field to 
a three-dimensional vehicle and obstacle field* This extension should 
include the testing of combinations of Obstacles which the vehicle could 
negotiate individually but possibly not in combination* From this will 
develop information which will define terrain conditions which pose to 
the vehicle true Obstacles for reasons of stability, power or torque re- 
quirements. An extension of this work will include soil-wheel traction 
parameters* This would define completely the terrain conditions and permit 
calculation of the vehicle ability to negotiate or not negotiate the in- 
dicated combinations of terrain geometry and conditions* 

A.l.d* Radio control - Mark Rinaldi, Mike Wolf 

Faculty advisor: Prof. 0* N. Sandor 

The objective of this task is to design and build remote control 
capability for the 0*4 scale model. This Involves establishing a radio 



link between the operator and the vehicle. The purpose of this remote 
capability is twofold. The first is that of observing the vehicle 
behavior under varying conditions for identical sets of instructions. 

The second is to provide a test bed for sensing and opera » t g systems. 
These objectives will ultimately require computer control of the vehicle 
through the radio link. 

The initial stage is the development of a manual input system. 
Additional features will be added, and refinements made until closed 
loop computer control is achieved. Incorporated in this task, and under- 
going parallel development is the vehicle function control. The received 
radio signal must be converted to a meaningful command to which the 
vehicle functions must respond. 

The parameters of the radio 1 ink system hive been defined and the 
system functional blocks designated, flexibility in the design of the 
system and its components has been stressed to accommodate additions and 
changes to the vehicle and the commands desired. 

The task of creating equipment designs for the functional blocks has 
begun. Designs for the manual ^nput device, command encoder, radio link 
and conmand decoder are complete. They have been built and tested, with 
the exception of the "clocks" requi ed at the transmitter and receiver to 
synchronize the data reception with the data transmission. 

Significant progress was made toward realizing only the first of 
the above objectives. After some thought it was decided that a digital 
communications system would be most appropriate, especially to simplify 
Interfacing problems with the computer. 

From the several types of digital modulation available. Frequency 
Shift Keying (FSK) was chosen, primarily for its high noise immunity and 
ease of implementation. FSK can only be used for low data rates, but 
since high speed was not a primary concern, a data rate of" 106 bits/sec. 
was chosen. 

Another factor in selecting a data rate of 106 bits/sec was selection 
of walkie-talkies to provide the necessary radio lick. While their sma' 1 
band width limits them to slow data speeds, their advanta b ss include their 
low cost, high flexibility, and the ability to use them in a two-way com- 
munications link* 

The time from September 1972 to January 1973 was spent in developing 
general design Objectives end system requirements. Decisions concerning 
signal type, modulation type, transmitter type and initial design con- 
siderations were pade during this time* 

Most of this effort was centered around the construction of the 
Cosaumd Encoder, a block diagram of which is shown in Fig. 10. The func- 
tion of the Command Encoder is to select data from either the manual 
controller or the digital computer and arrange it in a format that is 
suitable for radio t r ansmi s sion. The two major "blocks" in the Encoder 
are: 

1. CONTROL LOGIC, takas the input data and arranges it in a 
suitable digital format and generates the necessary timing 
signals for the transmission of the data. 
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2. DATA INTERFACE, converts the digital output of the 
Control Logic into the corresponding FSK tones. 

The Control Logic, Data Interface, manual control panel and the 
FSK detection circuitry which reconverts the FSK tones back into digital 
data have been designed and constructed. A more detailed discussion of 
these devices follows. 

Summary of Work 

The biggest problem encountered was the design of the Control Logic. 
The data format chosen consists of a 8-bit serial for each command trans- 
mitted. Since the computer provides data i.i this form, this choice was 
most convenient. The structure of the command work is as follows (see 
Fig. 11): 

1. The first bit, called the flag, is always a logical one 
and is used to mark the beginning of a word. 

2. Bits 2-7 are data bits comprising the command. 

3. The eighth bit is a parity bit. Even parity is main- 
tained and is used for error detection purposes. 

As shown in Fig. 11, the bit length is 9.4 ms. In addition, a 9.4 ms. 
timing signal called DATA SYNC precedes every command work and is used to 
synchronize the vehicle clock with the incoming data. A complete schematic 
of the control logic appears in Fig. 12. 

Definition of input symbols : 

1. DATA VALID - a high (logical 1) on this line indicates 
that the vehicle has accepted a transmitted command. 

2. M.R. - Master Reset - a low here resets all systems in 
t le Command Encoder. 

3. COMPUTER CONTROL - a high level indicates that the system 
is under computer control. All input data cot.es from the 
computer output. A low level enables manual control of 
the vehicle. 

4. MAN. XMIT (a&b) - the transmit command from the manual con- 
troller is entered here. 

5. CLOCK - a 106khz square wave derived from the 5.26 mhz master 
clock is used to clock data through the system. 

6. DATA IN - data is entered here from either the manual con- 
troller or the computer depending on the state of the 
COMPUTER CONTROL input line. 

Definition of output symbols: 

1. CLEAR - a zero level indicates system clear, (not to be 
confused with clear for individual counters or flip-flops) 
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2. XMIT - a one indicates that the system is transmitting 
and switches the walkie-talkie to the transmit mode. 

3. DATA SYNC - explained previously. 

4. HOLD - a one indicates that the systemis processing a 
command. In the computer control mode it will insl. .cl 
the compel >r to hold any new data until the system is 
cleared. 


5. ERROR - a zero indicates that there is a communications 
problem and the vehicle refused to accept commands. 

6. DATA OUT - eight 9.4 ms. bits that :omprise a command word. 

The control logic was built and tested and after some trouble- 
shooting it met its design specifications. Having completed the control 
logic, the next goal was the construction of the DATA INTERFACE. 

Three frequencies had to be selected to represent logical one, 
logical zero, and DATA SYNC. The frequencies were chosen to be V7ell 
within the pass band of the walkie-talkies and yet far enough apart 
to allow for efficient filtering. All three frequencies were derived 
from the 5.42 mhz master oscillator, (see Fig. 13). 

A block diagram for the DATA INTERFACE appears in Fig. 14. The 
digital waveforms from the control logic enables one of the three out- 
put lines (2048 hz, 2620 hz or 1310 hz) on which appear 0-5 volt square 
waves of the corresponding frequency. The lew pass filter is used to 
eliminate some of the higher harmonics. The high pass filter eliminates 
the d.c. component and, the new symmetrical signal, is fed to a buffer 
amplifier and finally to the walkie-talkie. The 4-volts peak to peak 
output signal is reduced to about 200 mv at the walkie-talkie. A com- 
plete schematic of the Data Interface is shown in Fig. 15 and typical 
waveforms are shown in Fig. 16. 

The next task that was undertaken was the construction of a COMMAND 
DECODER to reconvert the FSK tones back into a digital signal. The most 
important part of the Command Decoder is the tone detector circuitry, the 
block diagram of which appears in Fig. 17. A complete schematic appears 
in Fig. 18. The operation of the tone detector is relatively simple. 

An amplifier of gain 10 is used to boost the input signal (1C 1) 
and a high pass filter removes the dc component. IC 2 is a buffer 
amplifier used to provide a low source impedance for the filters. 

The active filters (IC's 3, 4, 5) ere of the multiple feedback 
bandpass type with Q ■ 10. Trim pots Rl, R2, R3 provide center frequency 
adjustment for each filter. 

Following each filter is a peak detector the output of which is 
proportional to the amplitude of the input sinusoid. 

IC's 7, 8, and 9 are voltage comparators used as threshold detectors. 
Trim pots RA, R5, R6 provide adjustable threshold levels for each of the 
detectors. The DATA and SYNC outputs are TTL and DTL compatible. 
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Summary of Radio Link 

With regard to the work completed so far, the major problem was 
receiver noise. To ensure proper operation of the vehicle, a low data 
error rate is essential. Although FSK provides good noise immunity, 
walkie-talkies are very noisy devices. Initial tests were disappointing, 
showing poor system noise immunity. After some adjustment of the thres- 
hold levels and better shielding techniques the performance improved 
significantly and is satisfactory. 

Major accomplishments during the past year were the completion of 
the CCMMAND ENCODER including the manual controller. The construction 
of the tone detector and subsequent testing showed that the system is 
capable of transmitting relatively error-free data. In addition, the 
present system incorporates enough flexibility to allow for many changes 
in system requirements. 

The problem of converting the received conroand into action by the 
vehicle was considered in conjunction with, but separate from the radio 
link. The system decided upon has a discrete command format, with the 
vehicle being self reliable (an open loop system). The received and de- 
coded command is stored and then converted to a speed reference in 
accordance with which the drive power is supplied. Figure 19 shows the 
functional block diagram of the motor control system. The system pro- 
vided three areas of concern: a time base, a motor speed reference, and 

a motor driver circuit. 

The design of the motor driver circuit settled on a class "D" modula- 
ted pulse-width amplifier. This was selected because negligible energy is 
dissipated in the output devices since they operate as common-emitter 
switches. 

Suninary - Control System 

Progre- s has been made on the control system in the areas of speed 
control, command storage and conversion to analog levels. The system 
was designed to be an all purpose one which could be duplicated for control 
of each of the vehicle's four motorized functions. To complete the system 
reference, sensors must be placed on the vehicle and their output incorpora- 
ted into the control system. The radio link and motor control system were 
linked together and used to control the states of a vehicle motor. 

The mated systems were able to produce three output speed levels in 
both forward and reverse directions. Control was done remotely from the 
manual control panel. The clocks for this were signal generators which 
will be replaced by individual receiver and transmitter clocks when they 
have been built end tested. 

The test showed the feasibility of the radio control system in the 
actual control of motor states. Hhlle much remains to be done before the 
system is complete, initial strides have been made. 

The time-base for the pulse generator, was drived from the master 
clock in the receiver decoder logic. This was divided down to approxi- 
mately 600 hx. This, was than fed to a one-shot multivibrator that reset 












36 , 


the ramp generator periodically. The voltage ramo was generated 
using a constant-current source transistor charging a capacitor. 

The voltage rises as a ramp since: 

V 

cap 

with I c h ar g e constant. Components of the ramp generator (Fig. 20) 

were picked so that the ramp could reach 5 volts before being reset by 
the discharge transistor. Two ri ups were made available by the opera- 
tional amplifiers: the first used as a buffer amplifier with a gain of 

ten, and the second as an inverting amplifier. The ramp outputs were 
sent to the comparators in the driver circuit, (Fig. 21). With a zero 
reference levjl the outputs were +12 for one and -12 for the other. 
Since the ramp did not reset to zero exactly, a dead-band exists which 
gives a stop level for the motors. 

Once power control was established, a reference level was needed, 
Fig. 22. The data was to be provided by the radio link, and a separate 
storage register was necessary to hold the commands. The only commands 
were speed magnitude and direction. Output from the J-K flip-flops 
were routed to the appropriate summing resistors via the two input and 
gates. The operational amplifier acted as a scaling adder in which the 
reference resistors provided several magnitudes of current to the 
amplifier inputs. The circuit was a simple digital -to-analog converter 
shown in the table below. 


Ex. BINARY 
0 0 0 
111 
110 
10 1 
0 0 1 
0 10 
0 11 

Storage 4 5 6 


ANALOG 
0.00 
► 5.00 
+ 3.34 
+ 1.67 

- 1.57 

- 3.34 

- 5.00 


if 


i 


charge 


dt 


This analog output was then sent to the comparators in the driver circuit. 
The pulsewidth was determined by the length of time the ramp stayed below 
the reference level (above the reference for the negative driver compare* 
tor). Fig. 23. A shift of reference voltage would therefore alter the 
width of the pulses and generate a change of speed. 

For the future, the completion of the first stage will be the installs* 
tion of the complete radio system for open-loop manual control. This will 
be exceeded first to closed- loop control, and then to computer operation. 
The on-board motor controller will be completed for all of the motorised 
vehicle functions. Once the control of the vehicle has been proven with 
basic command inputs (l.a., forward, reverse, turn thirty degrees), a 
function controller algorithm will be prepared to respond to complex 
commands (l.a., f'*rn from heading of 300 degrees to a heading of 200 de- 
grees). 
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.2. Collapsibility and Deployment - John Almstead 

Faculty Advisor: Prof. G. N. Sandor 

A solution for the problem of stowage of the RPI-MRV in the Viking 
capsule, and its subsequent deployment on Mars, has been found. The 
third generation design had as one of its criteria, to fit in the Viking 
capsule and to deploy with the greatest degree of simplicity and relia- 
bility possible. 

The vehicle is "folded" in such a way as to locate the center of 
the payload package on the centerline of the capsule, and to deploy 
by use of the front pivot, frame pivot and rear strut pivots as des- 
cribed in the following paragraphs and Figures 1 through 8. 

The front frame articulation pivot is to be used for collapsibility 
and deployment only. Once the vehicle has landed on the Martian surface 
and has been deployed to its operational configuration, the front frame 
articulation pivot will be fixed. The folded frame will most likely be 
deployed by some form of stored energy, such as springs at the pivot 
points. Collars placed around the overlapping frame sections will make 
the frame structure rigid. 

The rear strut pivot and lengthening will be accomplished by a 
four bar linkage. The linkage was developed through the use of the 
complex-number technique of linkage synthesis for three precision points 
and specified fixed pivot positions. A computer program was written 
based on this technique and used interactively on a graphic time sharing 
terminal to synthesize and analyze the four bar mechanism that can be 
used to deploy the rear wheels. Analysis of the linkage transmission 
angles insured easy and positive deployment which can be accomplished by 
a stored energy system such as a spring. Once the rear wheels have been 
deployed and locked with a collar or clamp arrangement, they will remain 
in the operational configuration for the remainder of the vehicle life. 

The front axle-frame connection has been designed to permit collapsi- 
bility and deployment, minimize weight, and provide for vehicle syiranetry. 
The front axle and wheels may be rotated at the front axle-frame connec- 
tion to bring the front wheels and axle above the frame. This maneuver 
is accomplished by the same motor used to steer the vehicle, and permits 
vehicle symmetry about the plane of the vehicle frame. This allows the 
vehicle to perform the "emergency maneuver". Without moving the vehicle 
body, the forward frame is brought over the vehicle to the rear, and vice- 
versa, clearing antennas, scanners and the like. Thus reversed, the 
vehicle can drive away from a condition from which it could not have backed 
away. 


These three mechanisms required for collapsibility and deployment 
are being designed, implemented and tested on both the 0.2 scale model 
and 0.4 scale sadel. 

Evaluation (.f the collapsibility deployment scheme is necessary 
with reference to the requirements for operation. In the near future, it 
is expected that requirements of the rear strut and of the forward frame 
for unfolding and locking the linkages in operating position will be in* 
vestigated. The collapsibility scheme will be evaluated relative to relia- 
bility and the operation of the RPI-MRV. The method of deployment will 
be fully implemented on the 0.4 scale model. 
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A. 3. a. Wheel Tester and Grouser Design - J. Ozimek and P. Morino 

Faculty Advisor: Prof. G. N. Sandor 

The purpose of this phase of the project is to determine the per- 
formance aspects of the wheels of the Martian roving vehicle. The project 
was divided into fo.r related segments discussed below. Basic equipment 
consisted of a 24' long plywood box 1' high and 1' wide. A test vehicle 
held the wheel and a drive motor arrangement. When the box was filled 
with sand, the wheel was driven on the sand, pulling the test vehicle 
which rode rails mounted on the box. 

Sand Preparation : It was desirable to simulate the soil conditions 

which the wheel would encounter over a year's span on Mars. This means 
that all possible combinations of soil conditions could exist, from fine 
loose sand to hard bedrock. From the traction standpoint, loose sand re- 
presents one of the extreme, conditions. For testing purposes, sand between 
two and sixty microns was obtained from a local creek bed. The sand was 
dried to as nearly zero percent moisture content as was feasible under room 
conditions. 

The experimental soil conditions could have been varied by the addition 
of moisture and by compaction, as well as the addition of various rock sizes. 
However, since there exist few published results against which experimental 
results could have been compared, the soil conditions were not varied. It 
was felt that dry sand representing the worst condition, would yield con- 
servative results. More significantly, the desired result of a comparison 
of the various wheel parameters, i.e. grouser design, is more readily dis- 
cemable if the surface conditions are maintained constant. 

Instrumentation : Initially a sophisticated arrangement was proposed 

in order to perform slip tests. After spending considerable time evalua- 
ting various techniques, a simpler scheme was settled upon without affect- 
ing the results to be gathered. 

To obtain information about the pulling ability of the wheel and the 
amount of power required by the wheel, a weight and pulley system was made. 
Heights, suspended out a window, were connected to the test vehicle with 
a cord and pulley arrangement. A wattmeter was used to measure the input 
power to the wheel drive motor. This was compared to the theoretical power 
required to raise the weights. 

Wheel rotation speed and cart translatory speed were also measured. 

Both the circumference of the test wheel axle and a cart wheel were seg- 
mented into conducting and non-conducting segments with tape strips. 

Brushes were installed on each. When electrical signals were sent through 
the brushes, the nonconducting portions produced a series of pulses which 
were printed on a strip chart recorder. The resulting plot was easily 
translated to determine the wheel's rotation speed and the cart's trans- 
latory speed. Relative slip was easily determined from these results. 

It was also desirable to measure the retarding force necessary to 
cause the wheel to slip. To determine this, the cart was affixed to the 
end of the test box via a linear spring. The spring deformation was then 
easily translated into the retarding force. 
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Grouser Design ; Few actual formulas or guidelines exist for the 
design of grousers; most design has been accomplished by use of empirical 
formulas or off the shelf designs. These formulas give only qualitative 
ideas of how grousers will perform under various conditions. 

The actual material used for the grousers was old tire rubber. Ini- 
tially the pattern of Figures A.3.a-1 and 2 was cut. The .3" thick lugs 
were glued to the rim surface, covering 90 percent of the surface. A 
second design removed every other grouser lug resulting in 45 percent of 
the rim surface being covered. The third design consisted of attaching 
the removed lugs on top of the remaining lugs. This resulted in the pattern 
of design #2 with a lug height of 0.6 inches. A fourth design used a lug 
of the pattern as shown in Figure 26, covering 60 percent of the rim sur- 
face. On recent lunar missions, this pattern proved to be effective in 
use with lunar roving vehicle wheel j. A fifth design merely doubled the 
fourth design lug height from 0.3" to 0.6". The lug to surface ratio was 
not altered. 

Testing Procedures : Two types of tests, slip test and power test, 

were performed for the five grouser designs: 

Slip Tests: A linear spring attached to the test cart 

provided the pulling force. For six variac settings (i.e. 
wheel rotation speeds), pulling force readings were taken for 
initial wheel slip, and total wheel stall. Two surfaces: l'> 
deep dry sand and 2' wood, were used. Power readings were 
also taken simultaneously with all tests. 

Power Tests : Weights suspended out the window were con- 

nected to the cart via the pulley arrangement. Runs from nine 
to thirteen feet for load weights of zero to two pounts were 
performed. Input power readings were taken in all cases. 

Results : Figure 27 summarizes slip tests on the deep dry sand 

surface. Types 1, 4, and 6 yield fairly similar results. Type 2, al- 
though showing poorer characteristics at low speeds, proves to be superior 
at speeds likely in the desired operating range (approximately 2 mph). As 
can be seen in Figure 28, type 2 again proves to be slightly superior on 
a hard surface. 

Stall slip results are summarized in Figure 29. Under equal 
loading types 1, 4, and 5 dig into the sand to a depth of one to two 
inches at stall slip. Types 2 and 3 penetrate to a depth of up to 4. or 
5 inches. Type 2 proved to be the best of the five designs tested. 

Figures 30 and 31 show cart speed characteristics versus power input 
for both a no vertical load and a two pound vertial load test. Since ex- 
perimental error may account for some of the differences, it appears that 
the different grouser designs vary little in efficiency. All designs 
exemplify linear speed-power characteristics. 

Conclusions: 1) Grouser thicknesses greater than 0.3 inches (for the 
6.85 inch radius wheel) will not perform well. Although type 2 generally 
proved to be the superior of the five designs, other basic designs should 
be tested. 


★Zero vertical load means that the wheel supports only its own weight and the 
weight of its support lever and drive. 
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Figure 25. RPI-MRV 0.4 Scale Toroidal All-Metal-Elastic Wheel on 
RPI Test Rig* with 0.3" Thick 45X "Hash" Grousers 
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2) The wheel shows excellent flotation. *>en carrying a 28 
pound load (comparable in scale to actual weight), the rim does not 
penetrate below the sand surface except when slip occurs. 

P roposed Future Work : Since this test facility produces only 

results which can be compared only on a relative basis with other 
wheels tested in the same facility, it will be necessary to test other 
proposed wheel designs to determine how the RPI wheel design ranks 
relative to other wheel concepts. Additional testing will be required 
and should include variation of soil parameters, by, for example, in- 
creasing the moisture content of the test soil and running the tests 
over a hard ° iface. The solid state tachometers will be perfected to 
make the output of the tests more easily reduceable. 

The ultimate objective is to optimize the RPI-MRV Toroidal Wheel/ 
Grouser Design. 


A.3.b. Wheel Analysis - W. Erableton 

Faculty Advisor: Prof. G. N. Sandor 

The task of developing reliable methods for analyzing and designing 
wheels for the Mars mission has been pursued for several years, and a 
number of computer programs have been developed. However, none of the 
programs correlate effectively the experimental data for the 0.184 scale, 

0.4 scale, or full-size models of the RPI torsional wheel. The objective 
of this task is to obtain an effective analysis method which can repre- 
sent accurately the behavior of the wheel. An effective method will 
provide a rational basis for the design of wheels for defined specifica- 
tions. 

Previous analytical methods considered a finite pin-jointed rim, 
neglecting the bending moments of the curved rim. Current analysis re- 
gards the rim as a continuous deformable member. Deflection character- 
istics for a concentrated point load art considered. 

Initially an Iterative procedure wat developed. The wheel was 
visualized as two separate configurations: (1) System "A", Fig. 32, 

with an inelastic rim and a series of radial linear springs (the hoop 
"spokes"), and (2) System "B", Fig. 33, with an elastic rim with the 
"sp Aes" replaced with radial forces. The basic iterative procedure de- 
veloped as follows: 

1. Apply an initial load to system "A". 

2. Determine the resultant radial spring forces from "A" 
and substitute for the radial forces of system "B". 

3. Determine deflection of "B" and proportionally deform 
the corresponding points of "A". 

4. Re-evaluate the resultant radial forces of "A". 

5. Repeat steps (2) through (4) until the vertical force 
components are "sufficiently small." 


Solid inelastic rim 
System of radial linear springs 
External force applied at center 
System of radial reaction forces 


SYSTEM "A" 


FIGURE >32*. System "A" 




Elastic rim 


System of forces 

System of deflections due to system of forces 


SYSTEM "B" 

FIGURE 33. System "B" 



This iterative procedure was programmed and deflections for various 
load forces were solved on a digital computer. Results were compared 
with experimentally determined results for the 0,4 and full scale 
wheels. As Figures 34 through 37 indicate, the analysis closely pre- 
dicts deflections for small loads. However, large deflection analysis 
is of little value. As a result, a new analytical technique was sought. 

The analytical model described above predicts the deflection of the 
rim to be symmetric abqut the horizontal center line of the wheel cross 
section. Experimental results showed that the lower wheel half deflects 
considerably more than the upper half for both concentrated and distri- 
buted loads. The symmetric deflection is a result of the analysis of 
system "A". Therefore, the current analysis utilizes system "B" with the 
linear springs of system "A". The forces applied to the elastic rim of 
system ,, B n are transformed into an initial reaction at the base of the rim 
and a system of forces which depend upon the system of deflections. The 
final expression relating the point load and the resultant rim deflection 
is as follows : 


n P F 

y ± ■ k l 1 <5(4^ - a j)yj + c~ “ a F^ 

J ^ 

where y^ = rim deflection to i*^ point 
n *■ number of spring spokes 
Py = reaction force due to initial boad 
a F ■ angle on the rim where Py is applied 
Pj = forces applied to the rim 

= angle on the rim where deflection y^ occurs 
c ** spring constant of spring spokes 
E R = Young's modulus for rim material 
I R ■ second 

r * radius of the rim 
K = (Cr 3 )/(irE R I R ) 


P j - CX J 


«(♦ 


.-a.) - . I {[cos K (♦.-a.)l/CK?-l) 2 > 
1 3 K-2 J 


This expression yields a system of linear equations of the form 
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This system was programed on the digital computer, and the results 
compared with the experimental large load deflections. Figure 37 in- 
dicates that this current math model will be a valuable tool in future 
wheel synthesis. 

The current analysis now must be expanded to include the applica- 
tion of distributed loads. This may be accomplished by segmenting the 
lower portion of the wheel, and considering the contributions individual- 
ly from each segment. If the segments are made sufficiently small, 
distributed loads may be approximated. 

Proposed Future Work ; The extension of this work will be in the 
investigation of modelling the wheel with a distributed load on the rim, 
calculation of the resulting footprint and of the pressure distribution 
on the footprint. Once this is completed, the design parameters for the 
wheel will have been determined. The tools for design of the wheel, in 
its present configuration, to meet specific requirements will then be 
available. 


Task B. System Analysis 

The objective of this task is to develop a framework within which 
design decisions involving conflicting requirements can be made, and to use this 
framework to identify desirable and non-desirable features of a roving vehicle 
design. Progress for the academic year included identification of optimal 
vehicle system designs, an analysis of the sensitivity of some of these designs, 
a start on a program to observe the operation of these designs in real time, and 
an analysis of the MRV computer as a queue server. 

Progress achieved is described in detail under the following subtasks: 
1) System Modeling and Design Optimization, 2) Sensitivity of Optimal Designs, 

3) Time Simulation of Vehicle Performance, and 4) Specification of an Integrated 
MRV Controller. Finally, a projection of activities for the next period is in- 
cluded. 

B.l. System Modeling and Design Optimization 

M. O' Callaghan, C. Pavarini, F. Stolen 

Faculty Advisor: Prof. E. J. Smith 

The method of investigating properties of Mars-roving vehicle 
designs is covered in past progress reports, and specifically in Ref. 3. 
Briefly, the method involves three major steps: 

1. Identification of static system models for 3 cases 

a. 4-wheeled, direct coammlcating rover 

b. 6-wheeled, direct coomunicating rover 

c. 4-wheeled rover, communicating via a Mars 
orblter relay link 

2. Establishing a system evaluation function in terms 

of the model parameters (in this case, the product of 
science time and straight-line distance roved), 

3. Nonlinear progressing optimisation with varying limits 
on vehicle velocity, slope tfireshholds, internal tem- 
perature, etc* 
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The optimization process for the design of a moderately high 
order system, such as the Mars roving vehicle (MRV), is a process 
of specifying the information required by the subsystem designers. 

This information consists of the operating requirements to be met 
by the subsystem, and all of the constraints under which the design- 
er must work. For example, the designer of the communications sub- 
system would be given the information on the maximum weight, maximum 
power available and the required data transmission rate for a pulse 
coded modulation subsystem. 

Previous work in the systems analysis area has determined a 
reasonably good model for the MRV, Ref. 3. This model was refined 
and used to obtain several optimal designs for differe.it operational 
restrictions. These operational considerations include such considera- 
tions as the frequency of science stops and safety considerations such 
as the constraints on speed and slope which the rover would traverse. 

It should be noted that the model used includes not only the modeling 
for the subsystems but also includes the design tradeoffs and inter- 
relationships between the subsystem. Although these interrelationships 
are mostly of the form of power and weight distributions, internal heat 
dissipation is also included in the considerations for the thermal con- 
trol subsystem. 

The major area in which the final model differs from the previous 
model, Ref. 3, is in the thermal control subsystem model. In particular, 
it was decided that heat pipes were inadequate to maintain the internal 
temperature at 300° + 5°K and hence a compressor type of active cooling 
system was included in the model. Another change in the modeling was in 
the manner in which the batter would be used. Previously it was 
thought that the batteries would .>e used as long as the MRV was moving. 

In the present model, the batteries are used only when the rover exceeds 
a slope such that the radio- Isotope thermoelectric generators (RTSs) can 
not provide enough power for the rover's motors and accessory equipment. 

Once the model was refined to include the changes noted above, 
the equations which constitute the model were rearranged to allow the 
design states to be represented as combinations of a few state variables 
and all other factors which go into a design to be represented as func- 
tions of the state variables. Once this rearrangement was done, it was 
relatively easy to put the problem into a form such that the techniques 
of mathematical programing could be used to optimize the design of the 
vehicle according to some design criterion. There was some rearrangement 
and elaboration of the chosen design criterion in order to put it into 
a form which was applicable to the chosen state variable form of the 
model. 

Since not all designs would be acceptable, a reasonable set of 
constraints had to be placed upon the design states and upon selected 
functions of the states. Most of the constraints were of the form of 
non-negativity of various elements of the design. An example of this 
type of a constraint is one that requires the weight of the communica- 
tions subsystem to be positive. There are also other constraints on the 
system which are of the form of "total system" constraints and an example 
of which is that the vehicle and all equipment be less than the launch 
weight. Once all of the proper constraints are placed upon the design 
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process, the modeling is complete enough for the mathematical pro- 
gramming techniques of optimization. 

The routine used for the optimization was Fiacco and McCormick's 
Sequential Unconstrained Minimization Technique (SUMT). This was ob- 
tained from the Hartford Graduate Center of RPI in the summer of 1971. 
After some reprogramming of the subroutines which constitute SUMT, it 
was compiled under Fortran H and installed on disk. The routine works 
using first and second .partial derivatives or, with some degradation 
of efficiency, the first derivatives only or no derivatives at all using 
numerical differencing methods to approximate derivatives in the latter 
two cases. 

During the past academic year, the system optimization process has 
been completed for all three cases, yielding 36 designs for differing 
assumptions. To best illustrate the type and uses of the results, an 
example will be given. 

For a 4-wheeled, direct communicating rover with a slope threshhold 
of 20°, velocity limit of 1.5 m/sec, internal temperature bonds of 
295-305° Kelvin, and antenna diameter limited to 2 meters, the design 
determined appears in Table 1. 

The complete set of designs will be covered in a Technical Report 
to be issued in September 73. 

No further work is projected along these lines in the immediate 
futu.e. In the long range, however, it is expected that the method 
would be applied to the RPI-MRV as its design is advanced. 


B.2. Sensitivity of Optimal Designs - C. Pavarini, M. Ryan 

Faculty Advisor: Prof. E. J. Smith 

Consider the case where an optimal design has been determined for 
a given mathematical model of a system, and a given evaluation function, 
as in Task B. 1. Suppose that one or more of the components of this de 
sign must be perturbed when the system is constructed. This perturba- 
tion (s) may be due to factors not considered in the system model: 

1. Although modeling equations are continuous, perhaps 
only certain discrete values are feasible (e.g., it 
is probably infeasible to design a communications 
system with a bit rate of 1.38762 K, but 1.5 K might 
be acceptable.) 

2. Perhaps another constraint is added, a posteriori (e.g., 
it may be decided to reduce the slope threshhold as a 
safety consideration) 

3. Perhaps some off-the-shelf hardware can be used, but 
its parameters do not coincide exactly with the optimal 
values. 


Then, it would be desirable to know how to* Change the other components 
of the optimal design to maintain the optimal property of the design, 
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Subsystem 

Weight (kg) 

Power (watts) 

science 

. 70,1 

26.0 

COMMUNICATIONS 

* 

139.7 

176.8 

POWER GENERATION 

105.9 

- 

THERMAL CONTROL 

29.9 

day=220.1,night=395 i 

NAVIGATION 

15,0 

6.0 

OBSTACLE AVOIDANCE 

5.0 

15.0 

COMPUTATION, DATA 

46.8 

90.0 

VEHICLE STRUCTURE 

197,5 



Other design parameters : 

COMMUNICATIONS DATA RATE 
BATTERY ENERGY 
RTG POWER OUTPUT 
INSULATION THICKNESS 
PROPULSION POWER 
VEHICLE WHEELBASE 
PATH-LENGTH RATIO 
DISTANCE BETWEEN SCIENCE STOPS 
NUMBER OF RECHARGES/DAY 


14,847 KBITS/SEC 
7.60 WATT-HRS 
545,4 WATTS 
0.028 m 
165.1 WATTS 
3.12 m 
2.18 
3.9 KM 
1.0 (average) 
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subject to what is, in effect, a new constraint. This new solution 
has been termed the "perturbed-optimal solution," and methods for the 
solution of this problem have been investigated over the past two re- 
port periods. During the past academic year, the mathematical deve lop- 
men; of the method of solution has been completed, and application to 
the MRV designs has begun. 


The mathematical analysis work included several parts. First, the 
existence of a unique, continuously differentiable trajectory, para- 
meterized by the initial perturbation, that solves the sensitivity problem 
was proved. Second, a first-order approximation to the general nonlinear 
solution was developed. Third, error-bounds on this approximate solution, 
and methods to improve the accuracy of the solution were obtained. 

The mathematical results to date will be presented briefly here. 
Consider the nonlinear programming (NLP) problem 


minimize f(x) 

subject to gi (x) ^ 0 i - l, . . . ,m 

hj(x) - 0 j * 1 y 

where x is an n-vector and f, and the gj and hj functions are nonlinear 
scalar functions of the elements of x. At the solution to NLP, x*, 
there exist two other vectors (an m-vector, u*, and a p-vector w*) such 
that 


h.(x*) - 0 

u i* 8 l( x *) * 0 

g 1 (* # ) £ 0 

«,* >. 


j " 1 > • • • > P 

i m 1| « • « y ,ID 

i * ly • « • y Itt 

i ■ 1 . • • > m 


and Voftx* u*,w* » Vf (x*) li* 70,-fx'* 1 ) 

«* i ° 


J a) * 7 hj M * 0 


where V i< the gradient operator with respect to x. 


A "perturbed-optimal solution" problem is generated from (NLP) by 
forcibly perturbing one of the components of x* and seeking the solution 
to this new problem. The solution to the perturbed-optimal solution 
is denoted x' and can be found froai 


x' ■ x* +£x 

where the Kth component of x'Cx'k, k fixed, 1 k n) is now fixed 
prior to solving this new probleri by writing 

-'it * *\ * 1 \ 

When 

e) the functions f , (|.), (h.) ere twice continuously 
differentiable 1 J 
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4r ^ 

b) the solution (x ,u ,w ) satisfies certain 
sufficiency conditions for (NLP) 

r7 * . k 

c) the Rectors V f(x ), Vg. (x ) foi i such that 

gj( x ) “ 0, "V h . (x*), j £ and e K (an n- 

vector of zeroes except for a * 1 * in the kth place) 
are linearly independent 

d) g^(x*) ■ 0 implies u^* > 0, 

a linear approximation to the perturbed-optimal solution can be found. 

That is 


' + <Tx a. 


* k £ * 
x + ° * k 


(where >| *is a variable vector except for a * 1 * in the kth place) 
where i\ K is solved from 


*) - G* 


G* T ® 


H 


*T 


O 

O 


du. 

du) 

d * r 

4s. 

dx K 


* o 


ohere dio^fu.*), G w * •• > H** VK P C. 


The conditions a through d above apply to the MRV system design 
problem and model. The results generated to date are values of the 
"perturbation coefficients", the vector in the equation immediately 
above, which indicate the sensitivity of the optimal design to the per- 
turbed parameter. 


The application to the MRV design problem will be shown^by ex- 
ample. Place the optimal design (unperturbed) in a vector x . Then, 
mathematically, the new (perturbed) design can be written 

x • x* + 

where is the perturbation coefficient vector, and £x. is the magni- 
tude of the perturbation in the k-th component. Table Ii shows results 
for an initial perturbation in science welgu«. An interpretation of 
these results indicates that the optimal way to Increase (decrease) the 
weight of the science package is to increase (decrease) the insolation 
thickness on the equipment package, thereby decreasing (Increasing) the 
thermal control power required. This change can be reflected into a 
decrease (Increase) in RTG power required and thus * corresponding savings 
(expenditure) in weight, which can be put toward the change in science 
weight. Note that the solution accounts for the weight change in in* 
creasing (decreasing) insulation thickness, thus taking into account the 
coupling In the system design. I? should be not*d that these results 
hold only for the case described in detail in B. 1, and thus not only 
are there many perturbed- optimal solutions problems for each design (a 
new solution for each value of k), but that each different design has 



65 . 


It 

\ e 

1 E 


TAB IE II 

THE PERTURBED- OPTIMAL SOLUTION FOR A PERTURBATION IN WEIGHT 

OF THE SCIENCE PACKAGE 


EXPERIMENTAL SCIENCE TIME 
TOTAL SCIENCE WEIGHT 
TOTAL SCIENCE TIME 
COMM. POWER 
COMM. ANTENNA SIZE 
COMM. DATA RATE 
RTG POWER 
VEHICLE VELOCITY 
SLOPE THRESHHOLD 
THERMAL CONTROL POWER 
BATTERY ENERGY 
..OVING TIME/DAY 
WHEELBASE 

EQUIPMENT PACKAGE AREA 
INSULATION THICKNESS 
HEATER POWER 
7o TERRAIN IMPASSABLE 
PATH-LENGTH RATIO 
SCIENCE STOPS/METER . 
RECHARGE TIME 
LARGE SLOPE FACTOR* 


35.75 

1.0 

35.75 

0.0 

0.0 

0.0 

-7.27 

0.77E-02 

0.81E-19 

-7.55 

0.252 

0.082 

0.75E-04 

-0.62E-02 

0.55e-03 

-7.27 

-0.15E-05 

-0.39E-05 

-0.16e-04 

0.67e-03 

0.78E-03 


* Large slope factor is the inverse of the percentage of time the 
vehicle spends on slipes the slope threshhold (the optimal value is 
51 . 3 ) 
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its own set of solutions, A limited number of cases have been solved 
as of this date. All results will appear in detail in a September 73 
technical report. 

No further work on this task is expected in the near future. In 
the long run, it may be applied to the RPI-MRV as its design is more 
completely developed. 


B.3. Time Simulation of Vehicle Performance - R. Gallerie, T. Mans 

Faculty Advisor: Prof. E. J. Smith 

As the system models derived are static, they consider many effects 
in their average or worst-case sense. This subtask is an attempt to 
investigate these effects in time (or equivalently, in their time-varying 
sense), by simulating the performance of the vehicle designs in real 
(scaled) time. Some of the effects to be considered in this manner in- 
clude slopes the vehicle will traverse, failures of the long-range 
obstacle avoidance subsystem, such alternative concepts as regenerative 
braking, etc. In addition, the performance of the designs under more 
comprehensive and varied earth command can be investigated. 

The progress for this academic year has been limited to the develop- 
ment of a general purpose Monte-Carlo simulation program to allow this 
type of simulation. This has been accomplished since this subtask was 
initiated in February 73. The program includes the system models derived 
in B.l. modified to run "in time", and is modular in its construction to 
allow easy expanding of its capabilities. 

This program simulates the operation of a Mars roving vehicle, with 
particular attention being given to power consumption. The system model 
is an adaptation of the system model used in Ref. 3. 

The model looks at the rover operation on an instantaneous level, 
where the model of Ref. 3 was an average type model. The program keeps 
a running account of time, power be i”g ufid, energy level of auxiliary 
batteries, and distance traveled. It also calculates the total time spent 
on science experiments, the total time spent stopping to recharge the 
batteries, the total time speut in direct communication, end the total 
ideal distance traversed. 

The program consists of the main program and four subroutines. The 
main program reads the inputs, calculates any design parameters not on 
the input list, initializes all variables, and prints out the results. It 
also creates a boulder distribution on the sample roving grid, which will 

be explained in detail later. The program finally directs the coded commands 

to the proper subroutines. Each subprogram is a distinct mode of operation 
of the vehicle. The first subroutine is the recharge subroutine. It is 
called whenever the battery level is lower than a specified level. This 

recharge subroutine is the only mode which is not started by direct conmand. 

The second subprogram is the science subroutine.. It is called by direct 
connand whenever the vehicle is to stop and take pictures. It is not used 
for communication dealing with transmission other than pictures because 
the power consumption is not as great. The final subprogram, the one of 
most interest, is the rove subroutine. In this subroutine, the actual 
roving is simulated. The vehicle roves about on a "sample" Mars terrain. 


► 

w 
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During its rove, the vehicle searches for obstacles, picks the "best" 
path, and travels the desired distance, meanwhile keeping an account 
of the power and energy consumption. The most critical power loads 
are used at this time. This subroutine is also called by direct command. 

The Mars "terrain" consists of a two dimensional grid, 100 meters 
by 100 meters, over which the vehicle travels back and forth. The 
distance between grid points represents one meter. The obstacles are 
actually dimensionless points, located on the grid by a two dimensional 
coordinate system which is supplied by input data. The obstacles are 
given dimension by the CLEAR input parameter. This parameter specifies 
the distance by which the vehicle must avoid the obstacle noint. Larger 
obstacles are created by placing two or more obstacle points adjacent 
to each other. Boulder obstacles and crater obstacles are not differ- 
entiated on the grid. The obstacle profile of the grid has an impassable 
area comparable to an average Mars terrain; however, it may be changed by 
input data cards. The profile considers any boulder over .5 meters high 
as being impassable. The profile of small rocks (under .5 meters high) 
and slopes is not time invariant. Although it changes with time, it is 
repeatedly determined by equations representing the actual slope and 
boulder distribution of Mars. 

The equations governing power requirements are adapted from Ref. 3. 

An equation calculating the power to rove over small boulders is a 
second order polynomial approximation for the data from a program de- 
veloped under Task A.l.c. This program calculates power requirement as 
a function of boulder size for a Mars roving vehicle. This program 
assumes that the power to traverse slopes and the power to travel over 
rocks can be superimposed. 

The program features a special recharge option which is controlled 
by the input parameter NOPT. If NOPT equals one, any surplus power pro- 
duced by the RTG's will be used to recharge the auxiliary batteries. If 
NOPT equals zero, this recharge option is not used. 

It is planned that continued effort on the time-simulation of 
vehicle performance be applied with the objective of gaining additional 
knowledge regarding design trade-offs, and the relationships between 
static system optimal designs and the performance of these designs in 
time. 

B.4. Specification of an Integrated MRV Controller 

C. Boroughs, F. Samuel, J. Zeitlin 

Faculty Advisor: Prof. E. J. Smith 

It is the objective of this subtask to describe in detail the 
principal tasks the computer must perform when it is roving on the surface 
of Mars and to estimate the performance requirements of the computer im- 
posed by these tasks. In addition, this group provides support for the 
remote control of the RPI-MRV under Task A.l.d. 

A statistical approach using queuing theory and random variables 
was chosen because the main object of the study was the average behavior 
of the computer-vehicle syrtem rather than a second-by- second study of 
the vehicle behavior. The roving mode was chosen as the object of study 
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because of the major computing requirements imposed by the obstacle 
avoidance and guidance functions as well as the strict time constraints 
due to the desire for efficient roving. 

Of course, at this stage in the development of the MRV there are 
many unanswered questions concerning the implementation of various 
functions in the roving mode. Thus, while some systems such as the laser 
rangefinder have been well defined, others, such as the path selection 
and obstacle avoidance algorithms are still in the early stages of develop- 
ment. Still, some estimates can be made concerning the amount of computer 
processing required. The goal is to give subsystem designers of such 
software-bound subsystems as obstacle avoidance and path selection a kind 
of bound on the computational capability available to them. 

The computer-vehicle system interaction was modeled by assuming the 
computer could be viewed as a queue server, servicing "customers" on two 
queues, the fixed queue (containing normal ongoing functions) and an 
interrupt queue. By queuing theory, it was determined how the timing of 
service to various tasks constrains the cycle time of the ccai r -cer. 

Tables I T I, IV, and V list operation timings, fixed queue tasks and timings, 
and interrupt queue tasks and timings, respectively. Results using these 
assumptions have not been fully evaluated as of this date, but indica- 
tions are that for this model, a cycle time under 50 nicroseconds should 
suffice in the roving mode for the cases listed in Table VI. Further 
analysis of results will be completed by August '73. 

Support has also been provided to the vehicle group responsible for 
Task A.l.d. in defining the conmunication system for remote open-loop 
control and ultimate closed- loop computer control. 

It is intended that both thrusts of this task will be continued in 
the future. Performance requirements on-board computer will be 

refined as more reliable subsystem re ements are obtained and analyzed 
through queuing theory. The interface ,.ith the RPI-MRV group will be 
maintained and particular attention will be given to the implementation 
of remote c mtrol with the ultimate objective of achieving autonomous 
roving closed- loop control with a remote mini-computer. 


Task C.l.a. Laser Rangefinder: Pulse Amplif icatior. - D. Palumbo 

Faculty Advisor: Prof. C. N. Shen 

There have beeu five requirements placed on a laser Rangefinder 
used for Martian Terrain Measurements. 

1. Range + .05 meters 

2. Angle + 10 arc seconds 

3. Fast scan rate, on the order of 1 m sec. 

4. Low weight 

5. Low power consumption 

The first two requirements are to assure meaningful calculation of terrain 
data when determining the in -path slope* To simplify the calculations, a 
complete scan of the' terrain should be made in about 1 m sec to eliml..ate 
error due to vehicle pitch and yaw, which leads to the third requirement. 
The last two requirements originate with the intended use on a space 
vehicle. 
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• t c computer cycle time - variable ; 

t_ subsystem access delay - variable 
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t . laser range-finder settling delay’- 1.1 mil. 
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t * analog sensor . settling time - 20 microsec. 
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t* / 

♦ * • 

• • « 

Register-Register Instructions It- . 

* ^ 

(Add Register, etc.) 

*■ > . . • • 

* < * - • • . 

Register-Storage Instructions * 2t„ 

(Load frcm storage, etc.) • 

• • . • . * * . 

* . ‘ ^ ' 

4 ' . • * • * • 

Multiply ‘ . • 6t rt 

. . . • * * • . * c 

.(floating point or integer) 

. • « • 

• • • '* •' ’ '• • : 

Divide . . .. ’• 12t c 

• . * •* 

(floating point or integer) * 

' .... 

• • * * . • 

• • • • • • • 

. • 

• Sine, Cosine Tangent 100t c 



TABLE IV 


• FIXED QUEUE TASKS AND TIMINGS. 

♦ , • 

» * * 

■ * • 

1. ENG Gather and encode engineering data 

" • • .time= 611t 2t • "V.“ 

; c a 

» ■ • . : 

. 2 . POWER Monitor subsystem power consumption 

Time= 970t„ + J2t + 320usec 
»•.» . •• . c . c . ' 

5. TEMP Monitor subsystem temperature 

* time= 582t. +• 32t + J20usec 

C ci J 

• 4. SCAN Update terrain array 

time* 3266t c + 1446t a + 522,500psec 

5* ANNAL . Analyse array and find path . 
time= 110 t 711t c 

6* GENER Generate path commands 

! * time= 107t /% + 4t • 

c a .... 


7. SCIENCE Perform science monitoring 
time= 72t c + 4t & + 120usec 


+ 


' TABUS' V • ' ' ' * " * ! - 

. 4 INTERRUPT QUEUE TASKS AND TIMINGS * 

. . t • 

Interrupt Function Exec Time Arrival Rate 

• • • " • 

1* FCOM Fet -jh earth command 119t.+2t 35-0 

C cl . * 


1 . 


2a.TWARN 

Thermal warning 
• « •* 

253 t c + 8 t a +20 

0.0017 

2b.BWARN 

• * • • 

• • 

Battery warning 

209t +10t +40 

• C . cl 

0.0083 

2c.PV/ARN 

Power warning 

255 V 8 V 20 - 

0.017 

■ * 



' • • 

3a.0ST0P 

Stop Tor obstacle 

535V 9t a 

0.035 

3b.SST0P 

• • 

Stop for slope . 

537t c+ llt a 

0.035 

Aa.NTIMER 

• 

Normal- request for 

522t c »9t a ’ 

xo 

• • 

path command 

i 



4b.STIMER 

Special request 

56,590t c +638t a + 

0.07 

* 

• 

• 

after a reverse ’ 
• 

• 

296,000 


5. DCOM 

* . • « 

* 

Down link ready 

82t^+t ft 
c & 

1.0 

• 

for data 


■ 





I 

. i.. 


TABLE Vi ' ‘ 

LIST OF TEST CASES FOR THE MRV COMPUTER MODEL 


.• Standard run - subsystem access time (t a ) 'is 10 cycles,. 

* • * •• computer v/uits during scan * . 


-2. Standard run - t is assumed 1 cycle - , computer still 
• * a • • • 

‘ ' vaits- during scan ■ 




3- Standard run - t = 10 cycles, computer is assumed to' 

a. 

. execute othor tasks without waiting during 
, . . .. scanning . 


-4. Standard run - t =1 cycle, computer does not. wait 

St 

; • during scanning . . . 


8. Runs 5 through 8 correspond to runs 1 through - 4 
* . * 
except that all task processing times were doubled. 

. * t . , 

That is, the number of processing cycles were 
.doubled, but the number of subsystem accesses 
and sensor delays remained unchanged. 
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Previously, only the "Phase Comparison" method offered accuracies 
meeting the first requirement. However, this technique uses a heterodyn- 
ing process which requires prolonged focus on a single point. Semi- 
conductor lasers cannot be pulsed over .2 microseconds leaving gas lasers 
such as Helium-Neon for use. The phase comparison method then may not 
meet fast scan assumptions and will be heavier and more power consuming 
than a method using a semi-conductor laser. 

The timed pulse method uses a GaAs laser and can be made to yield 
a .05 m accuracies if pulse expanding techniques are used. Pulse expan- 
sion is used on commercially available counters claiming .1 nano second 
resolution. Such a system has been designed and is proposed. 

The time pulse method is shown in Figure 38. The operation of the 
system is simpler than the phase comparison method, but, unfortunately, 
not as accurate. Initial investigations showed that .25 naosecond time 
resolution is needed when measuring the elapsed time of flight in order 
to achieve + 5 cm. accuracy, Ref. 4. The system shown in Figure 38 was 
suggested by W. Kuriger, Ref. 5, and is reported to be capable of + 10 
cm. accuracy under laboratory conditions. Some modifications can be 
made to improve Kuriger 's design. 

The range is determined by monitoring a ramp generator (most likely 
with an A/D converter) . The amplitude that the ramp reaches is propor- 
tional to the length of time it is on. The ramp begins on a signal from 
the pulse generator and stops on signal from the returning laser beam. 

Figure 39 contains the modified circuit diagiam. One important de- 
sign change concerns triggering the electronics with on board and return- 
ing laser beams. Lasers have an indefinite delay time between pulse 
application and onset of the lasing mode. Triggering with a split laser 
beam will eliminate the error due to this time lag. The processing be- 
gins on signal from the storage-averaging branch. A two ma. pulse 
triggers the power supply to drive tne laser. The time delay between on 
board and returning laser pulses is recorded by the detector and output 
as a pulse of width, T. This pulse width is expanded and timed. The 
count resulting from the timing is used as a memory address for read only 
memories (ROM's). The range is output from the ROM's and stored to 
await averaging. The pulse generator is then signaled for another range 
beginning the cycle over. The member of ranges to be averaged depends on 
how much time is allowed under the 1 msec, scan requirement. 

The pulse width, T, proportional to the beam time of flight is pro- 
duced by an exclusive or gate as shown in Figure 40. The beams are de- 
tected by photodiodes, amplified and applied to the gate. The on board 
beam will be stronger and require less amplification, thus A 2 will be 
greater than A}. The function of the gate as shown in the truth table is 
to produce the pulse length. Since the gate triggers at a specific voltage, 
5 volts, the pulse length will also depend on the slope of the leading edge 
of the input pulse. This won't disturb the operation of the on-board 
circuit, but will affect the timing of the reflected circuit since the 
shape of this pulse will depend on distance traveled and type of terrain 
encountered. This has been Investigated, Ref. 4, and it was shown that 


dt - T l /(S/N)/2 


MIRROR 



TIMED PULSE 






FIGUHE 39 

MODIFIED TIME PULSE 











where dt is the timing error, T the pulse length, and S/N the signal 
to noise ratio. Using nominal values, a maximum erre:* of 3.75 cm. is 
calculated. Since this is 757. of the allowable error, every effort 
must be made to minimize it. Although the reflected pulse shape won't 
be identical to the output pulse shape, it certainly will depend on it. 

The laser output pulse *'.hape must be as square and as powerful as 
possible. 

The problem of extracting centimeter accuracy from the timed pulse 
method may be approached as follows, assuming rounded numbers for ease 
of calculation. The range depth of field, R, will be from 3 to 30 meters 
and can be related through the speed of light, c, to the time field by 
Tj *• 2r/c. Accordingly T. must be between 20 and 200 n sec. If a range 
resolution of 03 meters is required, a time resolvtion of .2nsec. is 
necessary. As stated earlier, soon to be available counters will have 
resolutions of 10 nsec. (lOOMhz.). A minimum magnification of 50 x is re- 
quired for proper resolution. 

The combination Rb, T, and C constitute a ramp generator similar to 
Kuriger's design. The difference here is the method of output. A large 
resistor is placed across the capacitor with the result shown in Fig. 41. 
V^ n is the pulse ir.put from the detector and is width modulated depending 
on range. Using a capacitor of .01 microfarads and a resistor of 470K 
ohms, the transfer function of Figure 42 was obtained. The slope, or 
magnification, of the system is 1300 for the first part of the curve, and 
300 for the second part. Note both slopes are well over one hundred. The 
resistor, R, may be replaced by a high input impedance schmitt trigger 
shown in Figure 43. The function of the circuit is to provide triggering 
at predetermined voltages, producing a square wave output which is more 
suitable for timing. Figure 43 shows the input and output wave shapes 
along with the circuit's transfer function. The transfer function shows 
the hysteresis exhibited by the circuit, allowing precise wave shaping 
and timing. 

Figure 44 shows the expander with R' replaced by a transistor cir- 
cuit. The two transistors are driven by pulses V}, from the detection 
branch, and V 2 » a pulse generator triggered by charges the capa- 

citor as before, but now the charge is metered out at a constant current 
determined by the transistor. As before 


Vc l " I 1 T 1 /C 

If the stored charge IjTi Is metered out at current Ij, 


Vc l " X 2 X 2 /C 


Then simply 


The capacitor would supply a constant current until it's voltage dropped 
to the point where the transistor saturated. The discharge would then be 
exponential (see V ). Little time was available to be spent on this 
circuit, but it shows premise for development in the future. 










FIGURE 44 - PULSE EXPANDER - METERED 
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Up till now the concern has been producing a pulfe width dependent 
on range and long enough to allow proper resolution. It is clear from 
Figure 42 that the pulse expander's transfer function is not linear, and 
computation cf range from timing the pulse width would require extensive 
calculation. However, the pulse can be timed and hardwired to produce a 
memory address corresponding to the correct range. 

Since T. , the unexpanded pulse width, is related to the range by 
T^ » 2 R/c, the field o£ will be 

20^1^ 200 nsec 


and 


Tmax.-Tmin. = 180 nsec. 

If the range is to be known within + .03 meters, a differential element 
of range, dR = .06 meters, would correspond to a differential time of 


dT^ = 2dR/c = ,4nsec. 

The number of unique time values needed to define the field is then 

H = 180/d^ = 450 

If a portion of the pulse expander's transfer function, say with 
slope = 100, was blown ”p, the relation would be similar to Figure 45. 

If Tj = 60nsec., the differential element would be from 59.8 to 60.2 
nsec, as shown. This transfers to T 2 as mean time T 2 = 6 microsec. 
and differential element from 5.98 to 6.02 microsec. If a 10 nsec, clock 
is used, 4 pulses will fit into iT^. 

T 2 max. will be given by 100 x T.max or 20 microsec. With a clock 
period of 10 nsec., this transforms to a maximum count of 2000 at T 2 <nax. 
The number of binary digits necessary to handle this count is 12. For 
the case being considered here, with T 2 “ 6 microsec., the count will be 
600. Notice that any count from 597 to 604 is representative of the same 
c-apsed time, T^, and the corresponding range since all are within the 
differential time element dTj, see Figure 45. 

The memory address is produced from the count by hardwired logic 
gating, Figure 46. A count of 600 is shown in the register. This corres- 
ponds to an elapsed time of 60 nsec, or a ran of 9 meters. This range 
would correspond ideally to a memory address of 150. When the count is 
anywhere from 597 to 604, 150 should appear in the memory address register 
through similar gating. If the entire transfer function is calibrated in 
this manner, a unique memory address will be produced for any elapsed time 
T^ within the differential limits. 

The memory component is the least complicated of the system. It's 
purpose is to store 450 values of range lor processing. Since the maximum 
range has been assuned to be 30 meters or *3000 centimeters, a 12 bit 
memory word length will be necessary, the whole memory consisting of 450 
words. 
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As far as the averaging process is concerned, every logic gate has 
associated with it a delay time. This time is constant and if accounted 
for will not hamper operation. The typical delay time for a gate is 10 
nsec. The data access time from memory is 100 nsec. Notice that the 
maximum output from the pulse expander will be 20 microsec., and if 100 
gates plus memory have to be passed, the delay from pulse expander to 
range will be 21.1 microsec. A safe figure for the total figure would be 
about 25 microsec. Since a GaAs laser is to be pulsed for 100 nsec, 
(maximum) at a duty factor of .1%, the minimum period allowable for laser 
operation is 100 micro sec. This can be lowered easily to 50 microsec. 
by shortening laser 'on 1 time, but clearly the laser is the dragging foot. 
At a period of 50 microsec., 4 points will be averaged 10 times in 2msec. 
meeting pitch and yaw assumptions. 

The major features of a power supply for laser operation are the high 
current needed for onset of lasing (10 amps'* and the narrow pulse width 
necessary to avoid burn out. At first a supply was designed for use with 
the phase difference method. The current pulse must have a flat top if 
the laser output is to be modulated. Transistor 2N5262 was chosen for its 
speed and current capabilities. The final design is shown in Figure 47 
with 15 transistors in parallel to supply 45 amps to the laser. The cir- 
cuit was never constructed because 2N5262 could not be obtained. The 
circuit would have been ideal for experimentation allowing low current 
(non-lasing) pulses of long length for heterodyning experimentation and 
short, high current pulses for time pulse experiments. When it became 
obvious that 2N5262 would not arrive, an alternate design was chosen. 

Figure 48 shows a power supply and the current waveform it produces. 
The operation is to charge the storage capacitor through the transistors 
until the capacitor voltage reaches V cc . SCR GA201 is triggered to allow 
discharge through the laser diode. This pulser was built with the wave- 
form resulting shown in Figure 49. Not_, this is not the current waveform, 
but the signal waveform received at the photodiode. The rise time of the 
pulse (100 nsec.) is far from the current rise time given for the circuit 
(20 nsec), Ref. 6 and 7. 

Recalling the requirements placed on the rangefinder, they are: 

1. Range + .05 meters 

2. Angle + 10 arcsec. 

3. Fast scan rate 

4* Low weight 

5. Low power consumption 

The timed pulse method meets all requirements and is recoomended for 
further development in two main areas of study: the laser power supply 
with related detector electronics, and the modified pulse expander in- 
cluding counters and .registers. 

Supplying 40 amps to the laser becomes a problem, but can be solved 
by eliminating all unnecessary components and wire lengths in the laser's 
circuit. Figure 50 sauws a suggested power supply scheme to eliminate the 
3 XN914 's as an example of power supply improvement. The detection may be 
improved by adding a differentiating op-amp. Ref. 8, to eliminate received 
pulse rise time and the error Introduced by this rise time. 
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The modified pulse expander is not a unique design. Hewlett- 
Packard counters model numbers 5360 and 5379 boast .1 nsec, resolution 
and use pulse expanding techniques. They cost around $6,000.00, There- 
fore further research into the modified expander is recommended towards 
a possible linear transfer function to simplify processing. If the pulse 
expanding technique can be debugged and made practical, the timed pulse 
method, as described, is a laser rangefinder for terrain sensing measure- 
ments for Mars application. 

It is planned to continue the design, construction and evaluation 
of this laser rangefinder alternative during the oming year. 


Task C.l.b. Laser Scan Methods and Phase Comparison Rangefinder - G. Herb 

Faculty Advisor: Prof. C. N. Shen 

This task addresses two areas of study: (1) a laser scanning system, 
and (2) the range finder electronics and methods. 

If the terrain data is to be collected as the vehicle is moving, the 
rocking and rolling of the vehicle as it travels the uneven terrain can 
be a large source of error. To attempt to reduce this error, the laser 
scan methods considered here are of a block scanning nature. That is, the 
entire field of view of the rover is scanned out very rapidly so that the 
rock and roll error between points is negligible and all points In the 
field may be considered to have been obtained simultaneously. 

The criterion for such a system are: 

1. No or few moving parts 

2. An angle incremented s. *• 'i.e. point spacings 

vary with distance fr . . vehicle) 

3. Scan requirements - * • . i-; ;1 angle 40° - 60° 

i i~.glc '*0 

A rapid scanning laser range finder these requirements was 

shown to be a feasible system. A non— of . » laser scanner using 

either digital deflection optics or a nu-- . > *. semiconductor laser diodes 

was found to be the most promising method complete range finder scanning 
system was proposed around this type of rapid scan. The proposed system 
is highly accurate because it does not use %nolog signals to determine 
the angles at which a laser pulse is transmitted. All directions have bi- 
nary labels which correspond to angles which can be known to a high degree 
of accuracy. 

The problem of phase blur due to the incident spot from the leser 
beam being distorted due to the beam's shallow angle of incidence was 
found to cause negligible errors compared with the accuracies sired for 
a reasonable range of beam diameter* 

A d*rign for a phase comparison, nodulated laser range finder has 
been developed which appears to provide e very promising way of achieving 
the desired + 1 centimeter accuracy at mid-range (■ 15 to 20 neters>. 


A section of the design was breed-boerded and tested. Although the 
tests were in good a gre e men t with predicted results, additional work will 



91 . 


be necessary to complete the laboratory testing of the entire system. 

Details of the progress are provided in the following sections on 
scanning systems and rangefinder system. 

Scanning Systems 

One possible system, the dual vibrating mirror scan, was, invesri- 
gatcd but proved to be inappropriate for a rapid scan system. Here, two 
more systems will be introduced which attempt to alleviate the problems 
of the dual mirror approach. 

Kuriger, Ref. 5, has worked on and proposed a hybrid scan system 
using a vibrating mirror for a vertical scan and a line of laser diodes 
to accomplish the horizontal scan. His rangefinder system is not capable 
of producing the necessary accuracy, but his scanning system does have 
some advantages. Kuriger' s system is presented with some minor modifi- 
cations in Figure 5i. The "pulse transmitted detector" has been added to 
conform with the information needed in one of the rangefinding schemes 
currently under study. In this system both the transmitter and the de- 
tector use the same vertical scan mirror. The diodes provide the horizon- 
tal scan, while the receiver looks at the entire horizontal scan area via 
the cylindrical lens. The use of the same scan mirror assures synchroni- 
zation between transmitter and receiver, and the narrowing of the detector's 
field of view to only one horizontal row at a time increases the signal 
noise ratio. Further advantages of this system are its compactness and 
small number of components. But although the horizontal transmitted an^le 
is known exactly, the vertical angle determination relies on the same type 
of mirror position detection as does the dual mirror system and suffers the 
same problems. 

An alternative system which seems to avoid all the major scan problems 
is proposed. Figure 52 shows the proposed system using digital deflection 
optics, but a matrix of laser diodes could be substituted just as well. 

The detection system *>sed is si liar o Kuriger 's except the detector 
section is mechanically a separate entity from the transmitter. When the 
mirror drive starts its vertical sweep, the synch (start) circuit detects 
this and initiates the transmitter's scan. The transmitter's vertical 
scan is approximately synchronized to the mirror, but the detector is given 
a wide enough field of view to compensate for a slight mis-match. Any mis- 
match which does occur does not continue to accumulate since the transmitter 
and detector are re- synchronized at the beginning of each vertical scan. 

The advantages of this system are many-fold. First, the actual laser 
scanning is completely solid state. Also the transadtted angles are known 
exactly, with each pair of horizontal and vertical having its own binary 
code label. Finally, although a hybrid detector scheme is used, data point 
accuracy is extremely improved because we are i ot relying on the mirror 
position to give point location, only to nar±.--'» the detector's field of 
view to the approximate azimuthal scan line of the transmitter. Th.'a system 
then seems to have the potential for extrersly accurate transmitted angle 
information while *Jir.. aining a reasonable signal to noise ratio at the 
detector by narrowing its field of view. 



KORIGBR SCANNING SYSTEM 
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One problem with all the angular scan methods is the incident spot 
distortion at the terrain. Since a collimated beam is transmitted at an 
angle, the actual shape and size of the spot will be distorted when it 
strikes all surfaces other than those to which it is normal (see Figure 
33). This is most severe at the maximum ranges. As can be seen from 
the figure 


W = — - — for horizontal terrain 

Sin ^ (i.e. in the plane of the river) 

Therefore, at 30 meters (Q = 5.7°) a 2.5 cm beam diameter (= 1 inch) 
would produce a spot length of 25 cm. This effect worsens for negative 
in-path slopes and improves up to a point for positive in-path slopes. 


This distortion results in different rays within the beam traveling 
different total paths. Thus a ray on the top of the beam in Figure 53 
would travel a total round trip distance of 2( A-l / 2) or At more than a 
ray on the bottom of the beam. This distribution of path lengths within 
a single beam will cause a time of flight ambiguity at the rangefinder's 
detector. In the case of a modulated beam type rangefinder which is des- 
cribed more fully in the following section, the phase of the modulating 
frequency on the returning beam relative to the outgoing * 2 am is what 
supplies the range information. Spot distortion will cam.? a phase blur 
of a magnitude proportional to the amount of distortion. Referring to 
Figure 53 and taking the speed of light to be 3 x 10 s meters/second 


A^ = ,024f (w 2 - d 2 )* 

hi 

where A $ = peak to peak phase blur in degrees 
w = incident spot length in cm. 

d = beam diameter in cm. 

f ** modulating frequency in MHz 

But since 


6 « <At)f 360° 

111 

where $ *• relative receiver to transmitter phase of modulating 
frequency 

Ji t • beam ideal round trip transit time in microseconds 

The error introduced by the phase blur is 

,024 (w 2 - d 2 )* 

-^T “ 360 At 

and is Independent of f^,. Therefore, the phase blur error is independent 
of the beam's modulating frequency and can be fixed by other constraints. 

Table VII shows values of as a function of range in front of 




the vehicle, R, and beam diameter, d. 


R 

A 0/0 

d = .5 

TABLE 
in % 
d = 1 

VII 

d « 2 

3 

.0131 

.0262 

.0524 

5 

.0158 

.0317 

.0635 

10 

.0178 

.0355 

.0709 

20 

.0183 

.0366 

.0732 

30 

.0134 

.0368 

.0737 


This can be done since R is a function of Q and w is a function of d 
and Q . Note that the errors introduced by spot distortion on a flat 
terrain are negligibly small over a large range of d. 

A Rangefinder System 

A modulated laser - phase comparison rangefinder was proposed and 
appears in block diagram form as Figure 54. The benefit of this type 
of system is that since the returned signal will be a sine wave of a 
specific frequency in noise, it can be put through a very narrow band 
pass filter to increase its signal to noise ratio. 

Heterodyning techniques (frequency translating) are used as a 
method of increasing resolution. A constant phase angle is maintained 
in a heterodyning procedure. Therefore a 45° phase lag at 1 MHz corres- 
ponds to a lag of 125 nanoseconds. After heterodyning down to 10 KHz, 
the 45° phase lag is maintained, but a 45° lag at 10 KHz corresponds 
to 12.5 microseconds. Thus a frequency decrease of 100 times has given 
a time expansion of 100. The proposed system also averages several cycles 
of phase differences to further increase the range accuracy. A complete 
description of thij system can be found in Ref. 9. 

The phase detection section of the rangefinder was bread-boaraed in 
the laboratory to check its operation. Due to equipment problems the 
heterodyning section could not be simulated. Without the frequency trans- 
lating section, the phase comparator had to be checked with a lower fre- 
quency. The time expansion feature was also lost. Phase delays of 
sufficient magnitudes for the comparator to resolve thus had to be 
artificially produced. These phase differences therefore represented no 
actual range measurements, but were for test purposes only. 

Since this self contained system had no significant noise, the Schmitt 
triggers were set at the more easily implemented threshhold of zero with 
zero hysteresis. In the actual system, more realistic values of the thres- 
holds should be chosen. The most logical place to locate this threshhold 
voltage at first glance, is at the point of maximum slope of the signal, 
since this would reduce the sensitivity of the trigger to premature or late 
switching due to noise spikes close to the threshhold. However, for sinu- 
soids this point is at zero voltage which eliminates the hysteresis region 
and it is the hysteresis that reduces the sensitivity to premature OFF 



100 

MHz 






FIGURE 54 


r*< 


96 . 








ST0H2. 


*» " " T V ' " ■ t* -yi 









switching of the trigger after it has just crossed the ON threshhold 
and vice versa (see Figure 55). Thus the actual hysteresis, and 
therefore threshholds, should be set as a function of the most probable 
noise amplitude of the processed signal. 

As for the system under test, all sections worked as expected. The 
output wave closely resembled Figure 56 (with the exclusive OR envelope 
missing of course). Since NAND logic was used, however, the wave forms 
were inverted. With the SAMPLE command held high, the clocked phase 
pulses were grouped into the appropriate number of samples for averaging. 
Varying the phase angle of the input signals varied the phase pulse width, 
and therefore the number of pulses in the output count, by a proportional 
amount. 

It is planned to continue laboratory development of the phase com- 
parison alternative and to compare its feasibility with the pulse ampli- 
fication system described under Task C.l.a. 


Task C.2.a. Error Analysis and Stochastic Fit for Gradients - P. Burger 
Faculty Advisor: Prof. C. N. Shen 

Because of the round trip communication delay of between nine and 
twenty-five minutes, it is essential that a Mars roving vehicle be equipped 
with autonomous terrain modeling and path selection systems. The objec- 
tive of this task is to determine the effects of terrain sensing instru- 
mentation errors on the estimation of terrain features, specifically 
slopes and gradients, and to evaluate the feasibility of employing a 
stochastic fit for gradients. These results can assist in defining terrain 
sensor instrumentation specifications on the one hand and can guide the de- 
velopment of path selection systems on the other. 

The terrain sensor system being appraised, Figure 57, involves a range- 
finder which is aimed along azimuthal and elevation angles to a point on 
the terrain. It is assumed that a scanning system based on increment in 
the elevation and azimuthal angles, Figure 58, will cover the interest of 
region in front of the vehicle. It has been determined that such a scanning 
and rangefinding system can be used to estimate the local gradient of the 
terrain. However, because of the sensor errors in range and pointing 
angles, the estimated gradient will differ from the actual gradient depend- 
ing on the magnitude of the sensor errors. It has been determined that the 
standard deviation of the gradient will be in the range of one to eight de- 
grees for standard deviations in rangefinding of 0.01 to 0.1 meters and 
elevation and azimuthal standard errors of one minute. The standard devia- 
tion of the gradient is also dependent on the gradient of the terrain, on 
the distance of modeled plane from the vehicle and the data point spacing. 

In this analysis, it has also been assumed that a rapid scan is used so 
that errors due to the pitch and roll of the vehicle can be neglected. 

Complete details of the analysis and interpretation are presented in 
Ref. 16. 


Mathematical Analysis 

The quantities R, 0 , and ^ are measured with respect to the co- 
ordinate system h", a", and b", fixed to the vehicle, Figure 57. With 
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Figure 59. RELATIVE IN-PATH SLOPES 
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The body-bound axis rolls with th'* angle $ and pitches through an angle 
£ about a reference frame h, a, and b formed by the local vertical 
and an axis in a plane containing the heading and the local vertical. 

The* coordinate transformation, Ref. 10, is: 


h 

a 

b 


C(^> B (f) 


h" 

a" 

b" 


( 2 ) 


where C(0) and B(? ) are the roll and pitch transformation matrices 
respectively. ' 


Determination of Slopes and Gradient 

A number of measurement points within a certain small area of sur- 
face, say 0.5m by 0.5ra, can be used to determine a plane in space, which 
may be written as: 


h « axj^ + bx 2 + x 3 (3) 


where x^ is the cross path slope and x 2 is the in-path slope, which are 
the two constant parameters to be determined, Ref. 11 and 12. The corres- 
ponding gradient of the plane, which is defined as the gradient of the 
terrain in that small region on the planet's surface, is 

Sg - (x x 2 + x 2 2 )* (4) 

If the slope Is less than a predetermined criterion, It is considered to 
be safe for the vehicle to travel ahead. 

In order to locate the plane by a number of measurement points, one 
may rewrite Eq. (3) as 

h i " a i X l + b i x 2 + x 3 < 5 > 

where h.,a., and b. are found from R. , O , and A in Eqs. (1) and (2) 
for each ith point. Theoretically, three ^points determine a plane (1*1,2, 3). 
For greater accuracy however, more than three measurements (i*l,2,3, . . . ,n) 
probably n*4 or 6) are needed to determine the slopes. A complete picture of 
the terrain in front of the vehicle can be constructed by modeling numerous 
adjacent planes, each covering a small area of surface. 

Perturbation of the Variables 

If a^ • and b^ “ S, are assumed to be true in Eq. (5) a least 
square error estimate cen be performed which minimises 


( 6 ) 
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where is the actual measured height and is the corresponding 
height in the modeled plane. 

In matrix-vector notation, Eq. (5) is: 

T. = Ax (?) 

The least square estimate of the parameter x. Ref. 13, becomes 

x = (A T A)" 1 A T ii (8) 

Perturbation of the Variables 
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The least square estimate above assumes the matrix A to be com- 
pletely deterministic, Ref. 14. In reality, however, there is error 
involved in the determination of a, b, and h, due to the errors in 
our measured values of <^ , Jj , R, £ and Q . 

If the symbol £ denotes a perturbation, then Sh, &a, and&b in 
terms of £$>, £ § , gfc i . (h . and St? are: 


(9) 
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la 

» D(h\a' , ,b ,, f 4^) 
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8* 

+ c(4)o(^)g(r,0,6) 
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• Covariance Matrix of the Variables 
The covariance matrix of the variables can be defined, Ref. iS, 


M = « 


fih 

6o 

ISbJ 


[Sh Sa Sbj! 


L 


(10) 


where E denotes expected value. 

If $($ , i; , , 8(3 and \9 are not correlated, then with the aid 

of Eq. (9) 
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Proa Eqs. (10) and (11) the standard deviations of h, a, and b, 
can be cooputed In terns of those of 0, f , R. , p and 9 for 
each point. These are known quantities that depend upon the accuracy 
of the measuring devices. 


fru^rlane* Matrix of the Slopes 


Eq. (5) can be written as 




i.03 


t 


h - Ax (12a) 

With h * h + Sh, A * A + $A, and x « x + . ’>) 

and Eq. (7) the covariance matrix of the slopes is 

e[$x$x T ] = F^:[ShSh T ] - 2[6AxSh T 3 - E(Sh(S3x') T ] + E^SAx^SAx) ^ r T (15 

This equation utilizes the least square estimate of the parameter 
x from Eq. (8) and the elements of the covariance matrix of the vari- 
ables for each point, which are determined from Eq. (11). 


Variance of the Gradient 


If the symbol <5sg denotes the standard deviation of the gradient, 

Sg, then from Eq. (4) 

-Yk “V\ 

dSg = (xj+ x\) Xjdx^ + (x*+ x^dx^ 

The variance of the gradient is: 

_x Xx l - + 7 « 6 ) 

°*a “ <Tk ' + ><:>** + ** 

where CT^ . £(&)’, * E<^\ <T^ . *>&> 

The covariances CT** t and 0*,x* can ^ found from E< 3* 

The value of tfsg is a measure of the estimation of the gradient. 

In general, there is a 687. probability that the actual gradient lies 
between the values of Sg + tfsg and Sg - oftg end a 95% probability that 
the actual gradient lies between the values of Sg + 2 tfsg and Sg - 2 tfsg. 

Minimum Variance Estluste 
A generalized form of Eq. (7) is: 

v - Ax (17) 


1 


where v is the noise due to measurement of h and A in Eq. (5). If 
the expected values of Eq. (17) are 

E(h) - K , E(A) • X (18) 

then E(v) - v»h-Ax-0 (19) 

by virture of Eq. (12c). Subtracting Eq. (19) from Eq. (17), with 
the aid of tv. definitions of Eq. (12b) t one obtains: 

v ■ (h • 6 Ax (20) 


Thus the covariance of v b ec om es 

R . E(vv T ] . E(Uh - &Ax)Uh - SAX)*] 

. - EgAxfch 1 ] - E^h(jASc)^] ♦ eJ£aS(&Ax> T 3 C. 21 ) 
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which Is the same as those terms Inside the bracket in Eq. (17). 
The cost function J can be expressed as 
J « v T R -1 v » (h - AxjV 1 (h -Ax) 

« h T R -1 h - x T A T R" 1 h - h T R -1 Ax ♦ x t A T R“ 1 Ax ( 22 ) 


The minimum variance estimate of x can be obtained by taking the 
minimum of J with respect to x. 


Min-J 


22 ■ -2A T R“ 1 h + 2A T R" 1 Ax 
dX 


0 


( 23) 


Thus the minimum variance estimate x becomes 


/V 

X 


(A T K“ 1 A)" 1 A T R" 1 h 


24) 


where R can be obtained from Eq. (21). 

The minimum variance estimate given by Eq. (24) is a weighted 
least square estimate. It will provide a better estimate of the 
gradient than the standard least square estimate of Eq. (£) ’ 'use 
it gives correct weight to the individual measurement poi- 

However, the numerical results in the example compu .*om 
Eq. (24) is very close to the results given by Eq. (17) . »use the 
weighting factor in the covariance matrix R is nearly proportional 
to an identity matrix. That is to say, the weighting factors in the 
diagonal terms of R are equal. This is equivalent to an unweighted 
least square estimate. 


Numerical Results 


It is assused the rover will use a split beam where A/3 is the 
difference in the elvatlon angle between the two beams and Ad 1* 
the difference in azimuth angle between any two laser pulses as shown 
in Fig. 58. 

The values of £ and & for each data point, along with the magni- 
tude of the cross-path and in-path slopes of the terrain, determine the 
data point spacing and the distance of the data points from the vehicle. 
After a number of data points are measured, they are transformed to the 
non- rotating reference system by Eqs. (la to 2c). Than a least square 
estimate of the slopes is obtained by Eq. (8). This least squa s estimate 
depends directly upon the values of h, a, and b for each measurement 
point. Utilising the value of the vector £ found in Eq. (8), the values 
of $>, | , K, p , and © for each data point, and the standard de- 
viations of those quantities, the magnitude of the gradient and its 
variance can be calculated. 

The standard deviation in gradient, okg, can be determined by 
utilising Eqs. (11), (U), (15), and (16). The variance of the gradient 
depends upon the cross-path and in-path elopes of the terrain, the roll 
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and pitch angle of the vehicle, the data point spacing, the distance 
of the data points from the vehicle, and the accuracy of the measur- 
ing devices. 

If the rangefinder scan rate is on the order of milliseconds, 
which is perfectly feasible with electronic scanning, then each set 
of 4 adjacent data points is measured practically instantaneously 
since the rover motion is on the order of seconds and all 4 points 
will retain the same r'-’atr.ve position to each other when they are 
transformed from the vcniclf coordinate system to the fixed system. 
Consequently, the rover can model the planes in the vehicle coordinate 
system, and then transform the planes to the non-rotating frame. Since 
the maximum positive or negative slope that the vehicle can navigate is 
± 25?, relative slopes might be as high as + 50°and still be navigable. 
Fig. 59. 

Also due to rapid scan, the error introduced by setting « 1° 

(which is a reasonable estimate of these quantities due to the constant 
rock and roll of the vehicle as it tra ,,n rscs the surface) will be con- 
sistently between 1° and 2° because it ^nly involves the error in trans- 
forming the already modeled plane from the vehicle systeir of reference 
to the fixed system. Consequent l’*, the effects of Vp and are not 
considered in the analysis. 

If a flat terrain is assumed, » 1’ and <Tq. » 5 cm, and 

A& and (and consequently the data point spacing) are varied for 
each set of 4 points, the graphs shown in Fig. 60 are obtained. Each 
solid line represents constant values for Ad and A 3 . At and Ab. 
are the data point spacings along the in-path and cross path directions 
respectively. For any constant value of A 3 and AS, the data point 
spacing decreases very rapidly as the scan approaches the vehicle and 
consequently the error in gradient rises very rapidly. 

In choosing an 'optimum' data point space it must be noted that 
an increase in the spacing decreases the error in gradient. However, 
this renders the data less meaningful, as more terrain has been over- 
looked. At the 30 meter range, a general picture of the terrain with 
a spacing of 2-3 meters is sufficient. At close range, (4-7 meters) 
the points should be at least as close as 0.66 meters because this is 
the width of the widest navigable t revice. In Fig. 60, the dotted line 
represents dbg vs. distance from the vehicle for such a scheme, where 
the data point spacing varies from about 0.5 meters at a distance of 4 
meters from the vehicle to about 3 meters at a distance of 30 meters 
from the vehicle. 

Each of the plots in Figs. 61 and 62 utilize the same optimum 
scheme as the dotted line in Fig. 60. Once again flat terrain is 
assumed for all the data points. In Fig. 61, dfe is kept at 5 cm and 
and 0 ^ are set to 0, .017 (l 1 ), and .1 (6'). The error in gradient 
is markedly reduced by decreasing o'* and d$ from 6' to 1' but any 
further increases in accuracy will result in bulkier and heavier equipment 
which is not very beneficial, ty comparing the plots for • ofc ■ 1' 
and trjf - o£ ■ 0° it is seen that a standard deviation in djg and dg of 
1' has very little effect on the error in gradient. 


6 . 
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In Fig. 62, Or is varied from 1 cm to 10 cm while 0J3 and 
are kept constant at 1' (.017). The error in gradient increases 
rapidly at close range for each value of o^. 

Fig. 63 plots a' s „ vs relative in-path slopes of from -30° to 
+50° at 4 meters from the vehicle. In all cases, the cross-path slope 
of the measured plane is 0° and there is at least 1 data point with a 
height h of 0 meters in each group of 4 points. The graph only extends 
to -30° because the rangefinder beam cannot 'see' over a negative in- 
path slope greater than this value at 4 meters distance. With = <J£= 1', 

Or is varied from 1 cm to 5 cm. A data point spacing of .6 meters is 
assumed. Fig. 64 plots Ogg vs relative in-path slopes of from -6° to 
+50° at 20 meters from the vehicle. All the other conditions are the 
same as in Fig. 63 except that the data point spacing is 1.2 meters in- 
stead of 0.6 meters. Finally, Fig. 65 plots graphs of o^g vs relative 
cross-path slopes of from 0° to 50° for distances of 4 meters and 20 
meters from the vehicle. In this case, the in-path slopes are set to 0°. 

Once again = og =1' and is varied from 2 to 10 cm. The data 

point spacing at 20 meters from the vehicle ic 1.2 meters and the spacing 
at 4 meters from the vehicle is 0.6 meters. 

From Figs. 63 to 65, it can be noted that the higher the relative 
in-path or cross-path slope, and the closer the data points are to the 
vehicle, the larger the error in gradient. All other factors considered 
equal, the error in gradient is larger for a given in-path slope than 
for a cross path slope of .he same magnitude. It is obvious in Fig. 63 
that Og should be as close to 1 cm as possible in order that o^ will 
be within 2° at 4 meters distance and high values of relative in-pfth 
slopes. 

This analysis, which is summarized in Ref. 17, leads to 3 very im- 
portant conclusions. First, the very large error due to the motion of 
the vehicle can be eliminated by utilizing a rapid scan rangefinder. 

Secondly, if Ofe and og are reduced below 1 arcminute, no significant 
reduction in error is obtained. Finally, most of the error in gradient 
estimation is due to 0 ^. The vari. .ice cf the gradient should be as low 
as possible, but if an upper bound of 2° is set, then Or must be as low 
as 1 cm. 

Due to surface irregularity, data points will not fall at regular 
intervals as they would on a flat plane, but will be scattered about un- 
evenly. A procedure must be developed to determine how each group of 4 
points is to be chosen prior to the modeling procedure. Another question 
that should be investigated is whether the modeling procedures for 
terrain at the 20 - 30 meter range should he modified so that the rover 
can model hills and craters as well as flat planes. More generally the 
problem is one of determining a complete picture of the terrain from an 
array of different sized and shaped planes, each of whose gradients have 
a different magnitude and direction. 

Task C.2.b. Gradient Estimates from Stereo Measurements - W. J. Pfeifer 
Faculty Advisor: Prof. C. N. Shen 

The method of calculating a stochastic estimate of the gradient 
of the terrain in the path of the Martian vehicle, Task C.2.a. has been 
applied to two stereo measurement systems. One system, called stereo 
angles, uses three angular measurements in elevation and azimuth to 


locate a terrain point. Another system, called stereo range uses 
three range measurements. These two systems are compared with the 
non-stereo system. Task C.2.a. and Ref. 16, and with the limits of 
present laser rangefinder technology. 

The standard deviation in the four point estimate of the gradient 
depends on the distance of the four points from the vehicle, the posi- 
tions of the terrain points in relation to each other, the gradient of 
the terrain and the measurement errors. Vehicle motion effect on the 
standard deviation in gradient has been reduced to an additive factor 
as a result of the assumption of a rapid-scan laser. 

The effects of the above factors on the standard deviation gradient 
are shown in Figures 68-71 for a stereo angles measurement system which 
uses 3 angular measuremen:s and in Figures 72-74 for a stereo range sys- 
tem which uses 3 range measurements. A desired maximum allowable standard 
deviation in gradient has been set at 2°. These results are summarized 
and compared with the non-stereo system which uses one range and two 
angular measurements. Ref. 16 and 17, in Table VIII. 

Within the limits of present technology, the non-stereo system 
appears much more practical than the stereo range system. Laser range- 
finder technology has not yet reduced range standard deviation o^ below 
0.01 meter. In fact, a more reasonable expectation for range standard 
deviation may be from 0.02 to 0.05 m. The results for stereo range have 
been presented in anticipation of improvements in rangefinder accuracy. 

At longer range the non-stereo system s better than the stereo 
angles system. At close range, stereo angle, may be comparable to non- 
stereo if the standard deviation in the transmitter-receiver separation 
can be reduced to 0.001m. This conclusion is also based on the assump- 
tion that angle standard deviation will be about one minute of arc. 

Mathematical Analysis 

The basic method of estimating gradient from four terrain point 
measurements is outlined in Ref. 16. For this application, the follow- 
ing modifications were made. 

a) Stereo Angles 

The elevation angles CC and jS, and the azimuthal angle 
© are measured with respect to the coordinate sy: »em h", a", 
b" fixed to the vehicle as shown in Figure 66. A laser is loca- 
ted at point T which is 3 meters in height. The elevation angle 
of the transmitted light beam is measured at T fron the 
horizon to the terrain point U. The elevation angle oc of the 
received beam is measured at a point L meters from the top of 
the mast. Seme type of scanning detector is assumed in the 
measurement of the angle oc. The angle Q is the azimuthal 
angle of the transmitted beam. 

From the geometry, with a sensor mast height of 3 meters, 
the terrain point coordinates h", a"; b" are found in terms of 
oc , (2 , a, and L. 
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I Figure 66. Stereo Angles System Coordinates 
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The coordinates h", a" and b" are fixed to the vehicle and may be 
transformed to the h, a, b reference system as shown in Eq. (2), Task 
C.2.a. 


b) Stereo Range 

The quantities M, N and K are measured with respect to 
the coordinate system h", a", b" which is fixed to the vehicle 
as in Fig. 67. The measurement M is twice the range to the 
terrain point U from the transmitter and receiver at T on top 
of the sensor mast. 
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M = 2 R (2a) 

The measurement of N is the sum of the ranges from the trans- 
mitter to the terrain point U and from this point back to a 
receiver at a vertical distance L meters from the transmitter. 

N = R + P (2b) 

The measurement K is the sum of the ranges from the transmitter 
at T to the terrain point U and from this point back to another 
receiver at a horizontal distance J meters from the transmitter. 

K = R + Q (2c) 

From the geometry the coordinates of the terrain points 
are found in terms of M, N, K, L and J. 
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Figure 67. Stereo Range System Coordinates 
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Analytic Solutions for Estimation 

Let the symbol & denote a perturbation. If 5 9^, , 8e* , 5 ft 

, 80 and $L are not correlated the following relation is obtained. 
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The standard deviations of h, a and b can be computed in terms 
of those of 0, f , cx , & and L for each point. These arc 

known quantities which depend on the accuracy of the measuring devices. 


b) Stereo Range System 
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following is obtained by the same method 

used to derive equation 
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The standard deviations of h, a, and b can be computed In terms 
of those of §, M, N, K, L, and J for each point. Equations (4) 

and (5) Include the matrix G which involves all the measured quantities. 
Therefore, the standard deviations in h, a? and b depend on the loca- 
tions of the data points in relation to the vehicle as well as on the 
accuracies of these measurements. 

The slopes and gradient are calculated by the Eqs. (12) - (16) 
presented under Task C.2.a. 

numerical Results for Stereo Measurement Systems 
a) Stereo Angles Systems 

ft 

It la assumed that the scanning mechanism can provide 
any desired data point spacing at distances frees 3 to 30 
maters from the vehicle. Since the date point spacing also 
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depends on terrain irregularities, a flat terrain is 
assumed at the start. 

The distance L in Fig. 66 is assumed to be 1.0 meter. 
This is the most reasonable value for a 3 meter sensor mast 
height. A Larger value of L can giver smaller errors in 
measurement. But, as L is increased beyond one meter, the 
vehicle's ability to observe negative slopes is decreased. 
Also, if the receiver is lower on the mast there is more 
chance that the reflected laser beam will be blocked by 
other vehicle equipment. A rapid scan laser is assumed. 


Figure 68 shows a plot of standard deviation in gradi- 
ent Ogg vs distance, b", from the vehicle for flat terrain 
with 0° relative slope with 0^=0. 01 m. The dotted lines are 
for c^=0.005m. Clearly has an important effect on o^g 
for the stereo angles system, especially -^.t close range where 
decreasing o£ from 0.01m to 0.005m reduces the value of o^g 
by fifty percent. As expected, larger data point spacing gives 
smaller values for the standard deviation in gradient. At 
close range, data point spacing of 0.66m or less is required 
because 0.66 meters is the width of the widest navigable cre- 
vice. If <S^g is desired to be less than 2° it appears that 
must be about 0.005 meter or less. 


The standard deviation in the distance L ( 0^) may re- 
sult from changes in the sensor mast, expansion or contraction 
due to temperature changes, or from uncertainties caused by 
the laser equipment. 


Figure 69 shows the plot of versus relative in- 
path slope from -27° to +50° at a distance of 4 meters from 
the vehicle. The vehicle cannot measure negative in-path 
slopes of more than -27° because the terrain obstructs the 
line of sight. The cross-path slope, xj, is set co 0° and the 
in-path slope, X 2 « is varied from -27° to +50°. The standard 
deviations 0? and Cl are varied in Eq. (4) and the results 
used in Eqs. (15) and (16), Task C.2.a. to calculate the 
quantity o£ . In Figure 69 the data point spacing is 0.6 
meter as required at the 4 meter distance. The stereo angles 
system estimate of the standard deviation in gradient O^g 
is quite sensitive to the relative ln-path slope. 

Figure 70 shows the standard deviation in gradient o^g 
vs. relative in-path slope from -6° to +50° at a distance or 
20 meters from the vehicle. The terrain obstructs the view of 
the in-path slopes for more than -6° relative slopes at 20 
meters. For this graph the data point spacing is 1.2 meters. 

At 20 meters the standard deviation O^g Increases drastically 
with Increasing relative in-path slope. To obtain the quantity 
o£g less than 3° for all relative in- path slopes up to 50° 
the standard deviations 0? should be less than 0.001m and 
less than 6". 

Figure 71 shows the standard deviation in gradient vs. 
relative cross-path slopes from 0° to 4-50° for distances of 4 
and 20 meters with data point spaclngs of 0.6 meters and 1.2 
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meters respectively. If the standard deviations 0^*1' 
and sJ’l * 0.005m, or less, che value of O^g is less than 
2,5° st both distances. 

b) Stereo Range System 

As in the stereo angles system, it is assumed that 
the rapid scan transmitter at T, (Fig. 67), can provide the 
desired data point spacing for flat terrain. The separation 
distances L and J are both assumed to be equal to one meter. 

The one meter value is reasonable for the horizontal separation 
distance J because a larger J will require a more massive 
support. 

It is assumed that the standard deviations in L ard J 
are equal since the system configuration is somewhat symmetric 
in the horizontal and vertical directions. These uncertainties 
may result from changes in the support dimensions and from the 
laser transmitter and its 3 receivers. 

The range standard deviations ofo» o^, and cx R are all 
assumed equal to 2 o R . These standard deviations and the data 
point coordinates are employed to calculate the quantity c r 
in Equation (5), Eqs. (15) and (16), Task C.2.a. 

Figure 72 shows o' vs. distance, b*, from the vehicle 
for flat terrain with relative slopes » 0° and OU- 0.001m. 

The solid lines show the quantity o ' when the value of <Jj 
is 0.01m and the dotted lines when of ■ <rj is 0.005m. The 
standard deviations in gradient o' for the stereo range sys- 
tern are not sensitive to the standard deviations in separation 
distances o£ and cfj. 

Figure 73 is a plot of the standard deviation in gradient 
vs. relative in- path slope at a distance of 4 meters from the 
vehicle having data point spacing of 0.6m. In this case 
decreases as relative in-path slope increases to +50°. From 
the figure it 1 j concluded that to obtain the value of 
to be less than 2° for most relative in-path slopes, the magni- 
tude of must be 0.001m or less and the magnitude of o£ 
must be 0.005m or less. 

In Figure 74 the standard deviation in gradient is plotted 
against relative ln-path slopes at a distance of 20 meters with 
data point spacing of 1.2 meters. In this graph the value of 
the quantity o» increases slightly as relative ln-path slope 
Increases to +50®, To obtain a maximum standard deviation .°?g 
of 2° when the value of rt is 0.005m, the range standard devia- 
tion must be about 0.005m. 

The results which have been obtained are luaaarlxed in Table VIII. 
On the basis of these calculations and on the presented and projected 
state of the a.i; relating to rangefinder devices, it is concluded that 
the non-stereo range angle method studied under Task C.2.a. is superior 
to both the stereo-range and stereo-angles systems. 
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Task C.3. Path Selection System Evaluation 

S. Boheim, M. Szymanski, G. Ditlow, R. Campbell, R. bimonds 

Faculty Advisor: Prof. D. K. Frederick 

The objectives of the present effort have been threefo] the con- 
tinued development of a computer package which will provide the capabilitv of 
dynamically simulating a wide variety of sensor, terrain modelling, ar l path- 
selection algorithm combinations; the generation of criteria for the quantita- 
tive evaluation and comparison of the system's performance; and the use of 
these criteria in evaluating proposed systems. 

Previous efforts in the area of path-selection system evaluation have 
concentrated on the development of the major features of the computer simulation 
package and have been summarized in Reference 18. The progress attained during 
the past year consists of the following: 

a} Extension of the program to include an emergency backup 
mode and to incorporate nonzero vehicle dimensions. 

b) Addition of indices for safety and for obstacle detection 
errors to the quantitative performance measure. 

c) A comparison of the relative merits of vehicle- and 
vertical-fixed sensor configurations in terms of their 
effect upon path-selection system performance. 

In addition to the above tasks, considerable effort was devoted to the various 
and sundry chores associated with the enhancement and documentation of a large 
computer program. Such items included reprogramming for reduced running time 
.and storage requirements and improved output. 

After a brief description of the program's structure, the specific 
accomplishments noted above will be described, followed by a statement of the 
goals which have been set for the coming year. 

A. Program Structure 

As indicated in Figure 75, the computer program for the simula- 
tion and evaluation of proposed path-selection systems has been 
structured so that each block contains a separate function, thereby 
allowing substitution of alternate simulation schemes without major 
reprogramming. Briefly, the entire system operation can be divided 
into four major portions. These are: 1) terrain characterization, 

2) path- select ion system, 3) vehicle dynamics, and 4) system evalua- 
tion and simulation display. 

The terrain characterization portion is used to generate a very 
general terrain which includes boulders, crevases, craters, and cliffs. 

The path-selection system provides the functions of 1) sensor 
operation simulation, 2) terrain modelling (i.e., determining if an 
obstacle is present at a particular heading), and 3) path-selection 
algorithm. 

The third major portion of the program is the vehicle dynamics 
simulation In which the response of the vehicle to the commands of 
the path- select ion systemis simulated by determining how the moving 
Vehicle is affected by the terrain and how this in turn may affect 
the path- select ion system. 
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Finally, the system evaluation and simulation display supplies 
visual information indicating what the vehicle is doing, and what 
the terrain looks like. In addition, this portion assigns a numeri- 
cal value to the system's overall performance. 

In the following sections the major refinements to the computer 
package are described, followed by an evaluation of the effect of 
sensor gimballing upon the performance of the path- selection system. 

B. Vehicle Backup Capability 

In simulations reported in Ref. 18, no penalties were imposed 
upon a vehicle which struck an obstacle, thereby giving an unrealistic 
assessment of its performance. As a means of eliminating this defici- 
ency, an emergency path-selection algorithm with optional backup 
capability has been incorporated into the simulation package. This 
subset of the path-selection algorithm is only used when the vehicle 
encounters an obstacle it cannot negotiate. The vehicle backup 
function constitutes a separate block, Figure 76, which can be called 
only by the emergency mode of the path-selection algorithm. 

The two major considerations taken into account during the develop- 
meut of the vehicle backup capability were flexibility and realism. Any 
addition to the overall simulation program must not degrade these con- 
siderations. Flexibility is of primary importance because the substi- 
tution of alternative path-selection systems must be performed with a 
minimum of effort for the program to be useful. Realism is necessary 
to ensure reliable system evaluation. 

The backup block requires three parameters to be specified by the 
user, namely, the heading, speed (assumed to be constant), and distance 
of the backup maneuver. The program simulates backup motion and then 
passes the new vehicle location to the path- select ion algorithm after 
the backup is completed. Overall simulation data, such as total 
elapsed time, total distance traveled, and total battery time used are 
tabulated, thereby penalizing the quantitative performance evaluation 
figure computed by the program. 

The backup process is performed by the program in step increments 
of one half meter, unless otherwise requested by the user. The block 
also Includes the option of simulating an ideal mechanical sensor, at 
each increment. The mechanical sensor checks the safety of each incre- 
ment by calculating the in-path slope immediately behind the vehicle 
and comparing this slope with a prescribed limit (commonly chosen as 
30 degrees). If the calculated in-path slope exceeds the limit, the 
vehicle stops and does not traverse this increment. Control is now 
returned to the emergency mode for the vehicle to receive the next 
command. Figure 77-A shows the simulation of a safe backup path in which 
all slopes are within the allowed limit. In Figure 77-B the boulder is 
detected by the sensor since the angle O 3 is sufficiently large to ex- 
ceed the slope limit. 

% 

When the control returns from the backup procedure it is up to the 
emergency mode of the path- select ion algorithm to make the next decision 
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as to the vehicle's heading. Such a decision might call for turning 
the vehicle a certain number of degrees and scanning for a new path to 
the target or perhaps providing mission control a television picture 
of the situation and waiting for instructions. 

In the area of program refinements, effort has also been devoted 
to endowing the vehicle with a nonzero width so that it must have 
sufficient clearance on each side in order to pass by an obstacle. As 
originally programmed the vehicle was t.^ated as a point mass and in 
principle could pass between two very closely spaced boulders. 

A computer simulation demonstrating a vehicle backup maneuver is 
shown in Figure 78. 'When the vehicle fails to successfully negotiate 
the boulders a backup of six meters is performed. Upon completion of 
the backup the emergency mode commands the vehicle to turn 45 degrees 
clockwise and look for a safe path. This procedure allows the vehicle 
to find a safe path to the target. 


Performance Evaluation 

The performance index that had been used to date in the simulation 
consisted of an equally weighted sum of indices measuring: 

1) path length 

2 ) time taken to arrive at target position 

3) battery usage time 

The minimum values that could be obtained by a three dimensional 
straight line path from initial position to target were used as a 
base line for the performance evaluation. Subsequent simulations. 

Ref. 18, have shown that these measures alone are inadequate to evalu- 
ate the performance of a path selection system due to the fact that 
Increases in time and path length that occur when the vehicle strikes 
an obstacle are generally small and cause only small decreases in the 
performance index. Using these mersures alone near perfect scores were 
awarded to imperfect systems that failed to detect obstacles and ran 
into them. 

To improve the performance evaluation two additional indices have 
been added, namely 

4) safety index 

5) detection error index 

The first new index is an additional measure of path quality, 
penalising the system for approaching too close to and striking 
obstacles. The second penalizes the system for failing to detect any 
obstacles it may encounter in its path. The manner in which each index 
has been defined and implemented will be described below. 

1. Safety Index 

Having decided that a measure of path safety is desirable, 

the problem of measuring it for any arbitrary path on any arbi- 
trary terrain arises. The solution decided upon consists of 
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characterizing the location of each obstacle by a point in 
the horizontal (x,y) plane, computing the distance from the 
vehicle position to this location, then penalizing the system 
according to the distance between the vehicle position and 
the assigned obstacle location. The types of obstacles for 
which the safety index is applicable are boulders, craters, 
cliffs, and crevasses. For slopes which are too steep for 
the vehicle to negotiate no specific penalties are imposed 
by the safety index. 

To implement the safety index, the location of each 
obstacle ic cnaracterized by a point in the x,y plane, de- 
noted by the coordinates x 0 ,y 0 . For a boulder or crater, 
x 0 and y 0 are at the center of the obstacle's circular pro- 
file in the x,y plane. However, since the edges of cliffs 
and crevasses are second-order polynomials in the x,y plane, 
their location is characterized by that point on the edge 
which minimizes the distance in the x,y plane between the 
vehicle and the edge of the obstacle. 

The distance, denoted by D, in the x,y plane between 
(x 0 »y 0 ) and the vehicle's current position is computed for 
each obstacle at every sensor scan and when the ei-argency 
mode is activated. This distance is then converted to a 
danger cost according to the function sketched in Figure 78, 
where Dj is the radius of a crater or boulder, or half of 
the vehicle's length for a cliff or crevasse and D 2 is the 
desired pass buffer for the obstacle type. The values of 
the parameters Df and D 2 are specified by the user for each 
obstacle included in the terrain. 

Finally, a running total, P t , of the penalties assessed 
by the program is kept and the value of the safety index is 
computed according to the formula 

SAFETY INDEX - , * ■■ 


Clearly, the value of the safety index will be between 0 
and 1, with unity being the maximum attainable. 

2. Detection Error Index 

As indicated above, the purpose of this index is to 
penalize the path-selection system for failing to detect an 
obstacle which lies within its scan region during an« scan 
interval. Since the computer program knows the vehicle 
location and heading and the location of each obstacle in 
the x,y plane, it can compare this Information with the out- 
put of the sensor simulation to determine if a specific 
obstacle has been detected, l.e. , has caused the path- 
selection system to detect the existence of no-go lines of 
travel within the scan region. Failure to detect the 
presence of no-go lines of travel in the scan region when 
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an obstacle exists in the region is counted as a detection 
error and the system is penalized accordingly. 

The implementation of the detection error index is 
illustrated in Figure 79, where IND is a no-go line of 
travel indicator provided by the path selection system that 
has the value 1 when the system detects a no-go line of 
travel in the scan region and 0 otherwise. In the situation 
shown in Figure 79, the computer program has determined from 
the size of the scan region, the vehicle location and head- 
ing, and the obstacle location information, that a portion 
of the obstacle lies within the scan region during the 
current scan interval. Since the value o£ the indicator 
provided by the path selection system is 0 (indicating no 
detection of no-go lines of travel in the scan region during 
the present scan interval) the program has counted a detection 
error against the aystem. 

Finally, a running total, N DE , of the number of detec- 
tion errors assessed by the program is kept and the value of 
the detection error index is computer according to the formula 

DETECTION ERROR INDEX - fT^j 

DE 

Clearly, the value of the detection error index will be be- 
tween 0 and 1, with unity being the maximum attainable. 


The computer simulation shown in Figure 78 was run both with and 
without the implementation of the w-fety and detection error indices. 

The vehicle backup that occurs was initiated by an ideal mechanical 
sensor that detected the presence of a boulder in the vehicle path 
after the path selection system failed to do so on the previous sensor 
scan. Without the implementation of the new indices > the path selection 
system being used was awarded a performance index of 0.9/7. This nearly 
perfect score was awarded even thou s h tlu. system had failed to detect an 
obstacle in its path on the one occasion that an obstacle avoidance 
situation arose. 

The simulation was then duplicated, but this time the safety and 
detection error indices wen used in the system's performance evaluation. 
These indices were assigned weights of 0.70 and 0.15 respectively, while 
the path length, battery drain, and travel time indices were each assigned 
a weight of 0.05. A safety buffer of 1.5 meters from the edge of each 
obstacle was specified. The results of the simulation Indicate that the 
system was penalised for three detections errors and for approaching the 
boulders to less than 1.5 meters on three occasions. With the implement 
tat ion of the new indices, the same system performance was now awarded 
a performance index of 0.428. 

D. Sensor Orientation Comparison 
1. System Description 

One aspect of this year's work has been a comparison of vehicle- 
fixed and vertical-fixed sensor systems. In the former system, the 


sensor laser beam is maintained at a fixed angle of 82.4° with 
respect to the vehicle's mast whi.ch is itself fixed wi*h respect 
to the vehicle. Tn the latter case, the beam is meintained at an 
angle of 82.4° from the local planet vertical, see Figure 80. 

In both cases, 17 equally spaced range measurements are made 
within 40° scan ahead of the vehicle. For modelling purposes, the 
terrain is assumed to be linear tram t point underneath the vehicle 
to the impingement point of the laser beam. By comparing the re- 
sulting calculated slopes of the lines from the point under the 
vehicle to the 17 impingement points with acceptable maximum an** 
minimum slopes of + 30 degrees, the terrain modeller builds i 17 
element go/no-go i.iap corresponding to th** 17 possible i uths in ^ront 
of the vehicle. That direction closest to tne current target heading 
is selected for travel. 

If the slope is too great for travel in the direction the vehicle 
is currently heading, the surrounding paths are successively investi- 
gated until an unobstructed path is located and a steering command 
generated. Should all 17 paths indicate ^ostacles, i.e., slopes vrhich 
are too great, ‘■he vehicle rotates 30° clockwise and repeats its scan. 
The vehicle is initially pointed toward its target position and the 
mid-range sensor is used every 3 seconds. 

2. Terrain Description 

The two sensor systems were compared by usiog the computer 
program to simulate their responses over four terrains. In each 
case the terrain consisted of a sloping plane containing either a 
boulder or r crater directly between the vehicle's strrting p ;int 
and target. The terrain parameters are summarized in Table IX below. 


TABLE IX 

TERRAIN PARAMETERS 


Terrain 

Number 

In- Path 
Slope 
(degrees) 

Cross-Path 

Slope 

(degrees) 

Obstacle 

Type 

1 

- 25 

D 

Boulder 

2 

0 

25 

Boulder 

3 

♦ 25 

0 

Crater 

4 

♦ 25 

0 

Boulder 


3. System Performance 

Cases were first run in the absence of any random effects and 
then random noire was added to the deterministic values of the 
vehicle '<• in-path and cross-path slopes in order to simulate sensor 
attitude changes caused by passing over an Irregular terrain surface, 
e.g. , small roclu. and depressions. The terrain map and vehicle path 
ganerated by a typical run are shown in Figure 81, 
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Figure 81. Terrain Map With Vehicle Path 



While the terrain modelling subsystem and the path-selection 
algorithm used in these comparisons are not the most effective 
which can be designed and the terrains used are certainly not com- 
prehensive, certain conclusions can be drawn from the results. 

a) The vehicle-fixed sensor gave more reliable in- 
dications of positive obstacles than the vertical- 
fixed one. 

b) Both systems are sensitive to cross-path slopes, 
but the vertical-fixed version more so than the 
other. 

c) Negative obstacles, e.g., craters, appear to pose 
more serious detection problems than positive ob- 
stacles, e.g., boulders for both types of systems. 

d) Direct range measurements seem to provide a better 
indication of positive obstacles than the computer 
slope values used in the simulations. 

Summary 

Due to the nature of the computer simulation task which has 
been undertaken, it is necessary to continue to improve upon the 
program's efficiency and realism as it is being used to evaluate 
proposed path- select ion systems. Furthermore, as the program grows 
in size and complexity, the problems of documentation and handling 
grow. 


As more experience is gained in evaluating various path-selec- 
tion systems, the effectiveness of the present quantitative evalua- 
tion criteria must be scrutinized carefully and improved upon where 
deemed necessary. As a natural consequence of carrying out system 
evaluationr . efforts should be devoted to devising new path- 
selection systems which hopefully are capable of performing better 
than those studied to date. 

The activities of the coming year will be devoted to the tasks 
listed below, along with their target dates: 

1. Refinement of the simulation program: 

Addition of sideways roll (Dec. 1973) 
and wheel slippage (May 1974) 

2. Development of standardized test terrains (Dec. 1973) 

3. Study of vehicle behavior using only a short range 
(say 2 meters) senior (March 1974) 

4. Simulation and evaluation of the Cornell path- 
selection system (Dec. 1973) 

In addition to the above tasks, there are the ongoing tasks of im- 
proving the simulation program and evaluating alternative path-selection 
systems. 
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Task D. Chemical Analysis of Specimens 

One important phase of the initial missions to Mars is th« search 
for organic matter and living organisms on the martian surface. The present 
concept for attaining this objective consists of subjecting samples of the 
atmosphere and surface material to certain chemical and biochemical reactions 
and thereafter analyzing the products produced, probably in a combination gas 
chromatograph/mass spectrometer (GC/MS). The gas chromatograph is proposed 
for separating complex mixtures evolved from the experiments into small groups 
of similar chemical species. Chemical analysis of these groups would be 
accomplished in the mass spectrometer. It is the objective of this task to 
provide engineering techniques and criteria for designing such a system. 

The overall concept of the chemical analysis system is shown schema- 
tically in Figure 82. Most of the previous effort has involved an analysis 
of the gas chromatograph using simulation, Ref. 19 and 20. This technique 
uses mathematical models, which incorporate fundamental parameters evaluated 
from reported experiments, to explore various concepts and to direct further 
experimental research. The objectives of the task have been extended to 
include preliminary studies of the entire chemical analysis system. System 
topics being considered include carrier gas generation and removal, and 
limitations imposed by and upon the mass spectrometer. 

The task problems were attacked by a four-member team, each of whom 
pursued a specific assignment* 

1. GC/MS system concepts 

a. Mass spectrometer system characteristics 

b. Carrier gas generation and removal 

2. Chromatographic system analysis 

a. Multicomponent chromatography 

b. Chromatograph model improvement 

D.l. GC/MS System Concepts 

D, l.a. Mass Spectrometer System Characteristics - M. P. Badawy 

Faculty Advisor: Prof. P, K. Lashmet 

It is the function of the mass spectrometer to provide mass 
spectra for several inputs. Figure 82. This subtask had as its 
objective the characterization of the major properties of the mass 
spectrometer, the resulting relations to provide direction in de- 
signing an instrunent suitable for this space mission. In parti- 
cular, those characteristics influencing power, weight, and volume 
requirements of the system were of primary concern. Ref. 21. 

Prior work provided a parametric study of the coupling between 
the magnetic field and electrostatic potential in determining the 
ion masses to be detected. Instrument resolution and the nature of 
the ionizing process were investigated during the past year. 

Resolution, or the ability of the instrument to differentiate 
between similar masses at a given time, is estimated for the single 
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Figure 82 /. Schematic of Chemical Analysis System 







focusing spectrometer as follows, Ref. 22, 23: 


AM/M = 


s'+s" 


where 


+ 


Av 

V 


s',s" = slit widths at the source and detector ends, respectively 

a = magnetic field radius of curvature 
in 

V = voltage used for deflection of mass M 

V = voltage equivalent to range of masses emanating from 

the electrostatic field at different velocities 

M = dispersion in m .ss detected 


Similar but more complicated results for the double focusing mass spectro- 
meter have also been developed, Ref. 21. These results emphasize the 
importance of the slit widths, the radius of curvature of the magnetic 
spectrum analyzer, and the energy spread of the instrument. Table X 
summarizes the physical parameters of three mass spectrometers and pre- 
sents their respective mass resolutions as estimated from the above 
equations. Higher mass resolution can be achieved with more complicated 
designs, e.g., the toroidal electric sector field design. Ref. 24. 

However, for the current preliminary studies, the characteristics of the 
more simple, single or double focusing designs are satisfactory for 
representing the basic treatment. 

Instrument sensitivity or the minimum concentration level detect- 
able is determined to a great extent by the ionization process. The 
performance of the ionization chamber, b&sed on an electron impact 
source is characterized by the following equation. Ref. 25 and 26: 

I » b N O’ s I 

e 


where 


I, I e * ion and electron currents, respectively 
b * ion extraction efficiency 
M - molecular density 
& ■ ionization cross-section 
s ■ effective electron path length 

Difficulty in determining the ionization cross-section limited the utility 
of the equation. Ionization cross-sections for single electron emission 
ware available. Ref. 27. However, because knowledge of the components in 
the gas to be analyzed was uncertain, and determination of the ionization 
fragments was complicated, it appeared that theoretical treatment of the 
ionization process was beyond the scope of this preliminary study. 


A 



Experimental relations between the ion and electron currents will be 
sought to complete this phase of the task. 

Further work was concerned with determining the effect of major 
design variables upon power consumption of the mass spectrometer. 

Figure 83 shows schematically the relation between the various power- 
consuming elements of the instrument. Most recent studies involved 
estimates of power associated with the filament of the ionization 
chamber . 

In generating an electron emission current for ionizing the 
chemical specimens, the filament is a power consumer. Filament power, 
as determined by the filament current and its resistance, is aissi- 
pated through three mechanisms: 

1. Power associated with the electron emission current 
which is dependent upon filament temperature and volt- 
age. These effects are demonstrated in Figure 84 which 
is adapted from Ref. 28. 

2. Thermal radiation between the hot filament and the 
cooler envelope of the ionization chamber. This 
effect is proportional to the fourth power of the 
filament temperature. 

3. Molecular conduction by gas molecules in the space 
between the filament and the chamber envelope. This 
effect is proportional to the gas pressure. Ref. 29 
and, because of the low pressures of the mass spectro- 
meter, is generally of secondary importance. 

The steady state energy balance, obtained by equating the Joulean 
heating of the filament to the dissipation terms, was investigate! for 
several voltages and electron emitting areas. Table XI summarizes some 
of the studies and shows filament power is inefficiently used in pro- 
ducing an emission current. At the high filament temperatures con- 
sidered, most of the energy is dissipated as thermal radiation. Radia- 
tion shields may well improve this energy utilization. These data also 
show that fine filaments are required. 

These studies have determined the essential features of the filament. 
Future work will involve estimation of emission current ranges required in 
the mass spectrometer and the physical dimensions and voltage which re- 
sult in minimum input power. After these studies are completed, the 
accelerating chamber, the magnet, and the vacuum punping system will be 
considered. 

.b. Carrier Gas Generation and Removal - C. W. Jarva 

Faculty Advisor: Prof. P. K. Lashmet 

Carrier gas, which is relatively inert with respect to the samples 
being analyzed and which must be in supply for an appreciable part of 
the mission, is required to propel minute samples through the column 
of the gas chromatograph. On arriving at the mass spectrometer, which 
operates at a low pressure 10"^ torr), the major portion of this 



Figure S3. Power-Consuming Components of Mass Spectrometer 
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carrier gas must be separated from the sample and must be removed 
from the system. Thus the systems for generating, separating, and 
removing t.ie carrier gas are important to overall system design and 
operation. It is the objective of this subtask to generate the 
basic characteristics and limitations of several promising concepts 
to aid in the design of the final overall system. 

A study was undertaken, Ref. 30, to determine, as a function of 
the initial charge pressure, the physical characteristics of a 
spherical, high-pressure tank for storing the amount of carrier gas 
necessary for experimentation. Because of diffusion of the carrier 
gas, assumed to be hydrogen, through the vessel wall and its leakage 
through welds, valves, etc., Ref. 31 and 32, the gas storage capacity 
was greater than the requirement for the experimentation. Diffusion 
rates were determined from reported data. Ref. 33. To estimate leak- 
age through welds and valves, it was assumed the vessel would be 
subjected to helium leak detection when charged to a practical pres- 
sure taken as 34 atmospheres in this case. The minimum leak detectable 
by modem leak detectors, about 10’6 atm, cor /sec. Ref. 32, was extra- 
polated to account for hydrogen as the gas and for the actual charge 
pressure. The vessel wall thickness was estimated from standard 
equations, Ref. 34. 

Estimates of vessel size and weights as affected by charge pres- 
sure and operating time are shown in Figures 85 and 86. These estimates 
are based on the following assumptions : 

1. Construction - fiber reinforced stainless steel 
containing 1/3 metal by volume. Ref. 35. 

8 2 

2. Allowable working stress - 4.5 x 10 N/m 
(66,000 psi), Ref. 35. 

3 

3. Gas flow rate for experiments - 2 std cm /min, Ref. 36. 

4. Plight time to Mars - six montus 

5. Maximum temperature - 300°K 

Although diffusion and leakage were considered in these studies, these 
processes contributed only about IX to the vessel volume. The vessel 
weights, shown in Fig. 86 to be approximately independent of charge 
pressure, are for the shell only. The pressure regulator and its 
mounting boss will contribute an additional one kgm approximately. 

It was a major objective of this study to examine the merits of 
the gas storage vessel especially relative to a water electrolysing 
cell proposed earlier, Ref. 37. It is estimated that the gas storage 
tank will weigh approximately 60% more then a comparable electrolyzing 
cell system but will occupy about the sane apace depending upon the 
charge pressure. It thus appears that If reliable fiber- reinforced 
vessels having about four times the allowable working stresses of 
stainless steel, can be fabricated, the concept of the high-pressure 
storage vessel has promise, especially since the system Is simple. 



Figure 85. Effect of Charge Pressure and Operation Time on Size of Carrier Gas Sphere 
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Except for possible minor improvements in the calculational 
procedures, this preliminary study has been completed, and a technical 
report will be issued in the near future. 

A carrier gas removal system must be considered because the opera- 
ting pressure of the mass spectrometer is about 10 -7 times the pressure 
expected in the chromatograph. Ref. 38. The palladiun-silver alloy/ 
hydrogen gas separator has been proposed. Ref. 36. In such a system, 
hydrogen flowing in a palladium tube diffuses rapidly through the heated 
tubing wall in comparison to other gaseous materials. This subtask is 
concerned with developing technical procedures suitable for the design 
of such a ur't as well as with assessing its merits. 

For the initial studies, the steady state, of a system containing 
only hydrogen was considered, Ref. 30. Hydrogen diffusion rate through 
a palladium wall is proportional to the difference between the square 
roots of the hydrogen gas pressures at each face of the wall because of 
dissociation during diffusion. Ref. 39. From this diffusion rate equa- 
tion, a first order, nonlinear, differential equation was developed 
which related the rate of material flowing at any point in the tube to 
the rate the gas diffused out of the system. Combined with a second 
order nonlinear momentum balance equation, a describing set of equations 
for this system was formed which required numerical solution. This 
analysis assumed a flat velocity profile, pure hydrogen flowing through 
the palladium tube (no presence of sample), and radial flow of the 
diffusing gas through the palladium walls. 

These first studies showed the hydrogen pressure to be approximately 
constant over the tubing length, thereby eliminating the need for the 
momentum balance. The remaining rate equation was integrated directly. 
Figure 87 shows the estimated hydrogen flow rate as a fr<- tion of the 
inlet flow rate for such a palladium tube under the following conditions. 
Ref. 36: 

_3 

Inside radius of tube - 7.6 x 10 cm 

Thickness of tube wall - 7.6 x 10~ J cm 

Hydrogen pressure within tube - 4 atm 

Hydrogen pressure outside tube - 0 atm 

Temperature - 47 and 573 deg K. 

For this simplified system model, approximately all of the hydrogen 
would have diffused through a tube 0.27 meter in length. This compares 
to an experimental value of 0.31 meter. Ref. 36, an indication the 
model and diffusion coefficients reasonably represent the system. 

In an actual system, it is expected that the mass velocity of Figure 
87 will approach zero asymptotically because of the presence of samples 
which essentially do not diffuse through the tubing wall, m the future, 
equations accounting for the presence of sample will be developed. The 
present study also considered the hydrogen pressure outside the palla- 
dium tube to be approximately zero. Pumping means such as combining 
the diffused hydrogen with oxygen on the tubing wall. Ref. 36, is a 
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Flow Profile for Palladium Separator 



topic of future work. Other topics worthy of consideration include 
estimations of mechanical breakdown due to temperature and pressure 
cycling and chemical breakdowns due to contaminants. 


Chromatographic System Analysis 

a. Multicomponent Chromatography - A. J. Meisch 
Faculty Advisor: Prof. P. K. Lashmet 

Linear mathematical models of the chromatograph have been developed 
for single chemical components and have been solvea, analytically for 
an impulse input, Ref. 20, 40 and 41. Other work, Ref. 19, which used 
data from a chromatographic test facility, has shown that convolution of 
the actual input data with the theoretical normalized impulse response 
gives a reasonable representation of the actual chromatogram for a single 
component. 


It was the objective of this subtask to experimentally evaluate 
the mathematical model for multicomponent chemical systems using super- 
position of single component data. Specifically the equilibrium adsorp- 
tion model. Ref. 41, was considered: 


Y( 6) 

where 

0 

e 



^3Pe/4 7T 0^ exp £-(Pe/4£ 0) (0 

1 + iVmR Q ) 

dimensionless time - vt/L 
carrier gas velocity 
time 

column length 

the Peclet number, a dimensionless measure of 
sample diffusion in the carrier gas 

a thermodynamic parameter, specific to the 
chemical system considered 

composition response to a unit Impulse 


This model assumes the gas at each point to be in equilibrium with the 
solid adsorbent; i.e., the column is very long. The Peclet number is 
a function of the fluid mechanics and the physical properties of the 
system and is predictable. Ref. 20. The parameter mR^ is specific to 
the system and is determined from the system data using a curve fitting 
technique. Ref. 19. 


Extension of the model to multicomponent chemical systems using 
superposition of pure component data was proposed. Experimental work 
reported earlier. Kef. 20 on the binary system n-pentane/n-heptand on 
a mlcroporous adsorbent Chromosorb 102 coluom (see Table XII) showed 
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Table kll . . 
COLUMN CHARACTERISTICS 



' Chromosorb 102 .Car bow ax 1500 


length (cm) 

100 

100 . 

Outside Diameter (cm) 

0.32 

0.32 

Inside Diameter (cm) 

0.22 . 

% 

• 0.22 

Mesh 

•. 60/S0 . 

60/SO 

Particle Size (cm) 

0.025/0.0177 

0.025/0.0177 

Composition 

Microporous 
Styrene-Divinyl 
JBenzene Polymers 

Polyethylene Glycol 
20$ by weight 

Temperature Range 

Room-250°C 

Room-150°C 

Application 

» 

• ^ 

* • 

Separation of low 
.molecular weight, 
very polar, sub- 
stances such as 
water, alcohols, 
amines, acids, 
esters, ketones, 
ethers, aldehydes, 
and td a limited 
extent, light gases. 

Separation of high 
boiling polsT com- 
pounds . Gf eral 
purpose column for 
polar samples. 

* • 

• 

m • « 




superposition was only a first order approximation and large devia- 
tions were noted. This subtask had as its objective the verifica- 
tion of the previously observed discrepancies and the generation of 
additional multicomponent chromatographic data, Ref. 42. 

Chromatographic data for the system n-pentane/heptane on the 
porous Chromosorb 102 column have been obtained at 150, 175, and 
200°C and for mixtures containing zero to 1007. n-pentane by weight. 

As shown in Figure 88, the peak times for the multicomponent system 
lagged peak times predicted from superposition of the pure component 
data. This time lag occurred for both compounds for all mixtures con- 
sidered as shown typically in Table XIII and it increased as the com- 
ponent composition became more dilute. 

The peak time is primarily related to the parameter mR , Ref. 41. 
As the peak time varies with sample composition, it is possible that 
mRo is composition dependent. The component compositions in the helium 
carrier gas are not small (r- 10“^ mole fraction), and there may be 
possible interaction effects at the adsorbent surface. Figure 89 shows 
the values of mIU necessary to superimpose the model upon the data. 

This is empirical, and future work will consider interactions between 
the system differential equations. 

To determine if these interaction ef'ects were attributable to the 
porosity of the Chromosorb 102 adsorbent, the chromatographic separation 
of n-pentane and n-heptane was attempted on a nonporous Carbowax 1500 
column described in Table XII. Separation was not effected because of 
similar adsorption characteristics and large sample dispersion of the 
input sample. Possible modifications to the injection device to reduce 
the input pulse width as well as selection of a more favorable adsorbent 
are being considered. 


To improve data handling and consistency, a more automated data 
processing scheme outlined in Figure 90 vas considered. Implementation 
of the proposal with existing major hardware appeared feasible. How- 
ever, manual analysis of the stripchart chromatograms yielded about 
twice the resolution obtainable from the automatic processing system 
which was limited by the A/D converter of the available minicomputer. 
The proposal was rejected, and possible changes to the detector control 
circuits to obtain full scale readings with most samples are being 
considered. 

Experimental studies on the porosity and pore-size distribution 
of the adsorbents were begun. These data are required for accurate 
prediction by the improved mathematical models being developed (see 
Task D.2.b. ) . Initial data were obtained with a mercury porosimeter 
using a sample holder more suited to coarse adsorbent particles than 
to the 60/80 mesh particles of the chromatographic columns. The data 
appeared qualitatively reasonable, and work will continue upon receipt 
of a more suitable sample holder. 
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I • ‘ Flow Rate of 29.4 cc/min at 150°C 
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Composition 
weight ^ 

Time of Pentane 

Peak, sec 
* 9 

m 

Time of Heptane 
Peak y sec 

j 

100 H 

t - * 

• 

2*5.5 • ‘ 

I 

. 99 H - 1 P 

• 95.0 

29*1.0 

L* 

• 75 H - 25 P 

88.5 

307.5 

j r • 

.50 H - 50 P 

86.0 

327-5 


25 H- 75 P 

84.0 

• 30t.o 

i 1 

1 H - 99 P 

. 81.5 

432.5 

i 

100 P 

$ 1 
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Figure 89. Infect of Composition on Heptane mR ( 
for Pentane-Heptane System on C-102 
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2.b. Chromatograph Model Improvement - P. T. Woodrow 

Faculty Adviser: Prof. P. K. Lashmet 

Prior to July 1972 the mathematical model of the chromatographic 
column was based on a packing material which either was non-porous or 
allowed pore diffusion and subsequent adsorption to occur at such a 
high rate that the resultant system dynamics were not significantly 
affected. Thus, in modeling the system, only interp article and 
adsorbed phase mass balances were considered. The applicable partial 
differential equations, assumptions, initial conditions, and boundary 
conditions have been presented elsewhere. Ret. 2U. These equations 
showed the predicted chromatogram to be a function of three variables: 

mKg * a thermodynamic parameter, specific to the 
chemical system considered 

*= the number of transfer units, a dimensionless 

measure of diffusion rates to the adsorbent surface 

Pe ■ the Peclet number, a dimensionless measure of sample 
diffusion in the carrier gas. 

Solutions to this system of equations have been reported, and cases 
of equilibrium adsorption (Ntoc - * and non-equilibrium adsorption 
(N t oc finite) have been compared with experimental data, Ref. 20. It 
was concluded that for the ranges of parameters encountered in typical 
chromatographs, consideration of a finite rate of mass transfer (N tQG 
is finite) yields minor, at best, improvement in output chromatogram 
predictions. In addition, the predictions of peak height and spreading 
were inadequate for some exper cental conditions and chemical systems. 

It appeared that another mechanism not considered in the model formula- 
tion was responsible for the discrepancies between prediction and experi- 
ment. As the adsorbent (Chromosorb 102) was porous. Table XII, the 
mechanism of intraparticle diffusion warranted investigation. 

Consideration of the intraparticle diffusion process as well as 
the rate of surface adsorption yielded a new system of equations shown 
in Figure 91 and having the following parameters and variables not 
defined previously: 

m » thermodynamic equilibriun constant 

N Rn » the number of reactor units, a dimensionless 

RU measure of the rate of adsorption 

Pe. - Intraparticle Peclet number, a dimensionless 

measure of diffusion rates within the particle 

Pe_ » interparticle Peclet number, a dimensionless 

E measure of diffusion rates in the carrier gas 

L/R * ratio of coluan length to particle radius 

R_ » moles of fluid in particle .per mole of adsorption 

1 sites 

r - intraparticle space variable, dimensionless 
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Figure 91. System Equations for Inter- Intraparticle Model 


x * adsorbed phase co. >osition, dimensionless 

y * interparticle gas phase composition, dimensionless 

■ intraparticle gas phase composition, dimensionless 

y* * equilibrium intraparticle gas phase composition, 

dimensionless 

z ■ longitudinal position in column, dimensionless 

As an aid to envisioning the various transport processes incor- 
porated in this new chromatograph system formulation. Figure 92 is 
included. In the interp article region, material is transported by 
convection (i.e., by carrier gas flow) and by turbulent and molecular 
diffusion in the axial direction, which is characterized by the inter- 
particle Peclet group, Peg. The transport from the interparticle 
region to the intraparticle region is characterized by the dimension- 
less mass transfer coefficient, N t QQ. Transport by diffusion within 
the particle - the intrap article region - is characterized by the 
Intrap article Peclet group, Pe^. The rate of adsorption from the 
intraparticle gas phase to the adsorbed phase within the particle is 
characterized by N^j, the number of reactor units. 

The transfer function of the system of Figure 91 is given in 
Figure 93. Before attempting an inversion of this equation, a moment 
analysis of the model hereafter designated as the Inter- Intraparticle 
Adsorption Model was undertaken to determine gross behavior. 

The three moments of a chromatogram. Ref. 41, 



which arc derivable directly from the system transfer function end the 
moments of the input or sample injection pulse, give an indication of 
the time for peek appearance, the dispersion, end the skew respectively 



Ficure 92. 
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Tables XIV and XV compare chromatograph moments obtained from 
ethylene and acetone data, Ref. 20, with predictions from the inter- 
particle equilibrium adsorption model, Ref. 41 (see section D.2.a. 
above) and the new model. A priori calculation of (L/R), Pe^, Peg, 

Nf-OG was possible from thj experimental conditions. The values of 
itiRq determined earlier, Ref. 20, were used, and computations wl^ 
made for several values of Nrjj. In both cases the new model more 
closely predicts the mean /jj ^ although differences between the two 
models and the data are small. More importantly, the new modjsl is 
capable, in both cases, of better predicting the dispersion > 2 ’ and 
more closely approximates the skew JZ 3 . The values of N„tj necessary 
for matching jH 2 for both systems are of the same order of magnitude 
and are in reasonable accord with other independent data. Ref. 43. 

The differences in the gross behaviors of the two systems as repre- 
sented by the respective moments is primarily a result of the parameter, 
.tiRq as all other parameters are of the same order of magnitude. As a 
result of the moment analyses and parametric studies, it appears that 
the proposed model will yield more adequate chromatogram prediction. 

The solution to the equations of Figure 91 is presently being 
considered. It appears that inversion of the transfer function of 
Figure 93 is not feasible, and an analytical solution is probably not 
realizable. The simplified situation of negligible mass transfer 
resistance ) and equilibrium adsorption was also 

investigated. It appeared that the resulting transfer function while 
considerably less complicated, is also not easily inverted. Thus, the 
application of an appropriate numerical technique is required for the 
solution of the system equations. 

Prior to obtain the final numerical solution, several preliminary 
studies were required. The second order equations of Figure 91 require 
two boundary conditions which prescribe the solution behavior at two 
points in the space domain, as wall as a set of initial conditions. 

Thus a complete description of the interparticle mass balance of Figure 
91 becomes 


(l/Pe) ^y/«)z^) - (c)y 

7 <3z) - R = Oy/ ^©) 

( 1 ) 

y (z, 0 ) = 0 


( 2 ) 

y (0, ©) = A$ 

© > 0 

(3) 

lim y(z, ©) = finite 

e > 0 

W 


Z->QD 


where R is the rate of adsorption of the samples. Equations 3 and 4 
are the boundary conditions that have been previously used for develop- 
ment of the analytical models. However, use of Equation 4 as a terminal 
boundary condition is Impractical for numerical treatment an! must be 
replaced with an appropriate finite terminal boundary condition. It is 
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TABLE XIV 


MOMENT ANALYSIS AND PARAMETRIC 
STUDY - ETHYLENE 50°C. 


mR, 


0.194 


^1 , observed 

26.475 


predicted 

25.986 


( 2 ) 

**1 , predicted 

23.719 ‘ 



^2 

observed 


“(1) 

m 2 

predicted 


“( 2 ) 

M2 

predicted 


100 

.200 

400 

800 

1600 

3200 

6400 

12800 

25600 


7.024 


JL 


13.283 

6.973 

3.817 

2.240 

1.451 

1.056 

0.859 

0.760 

0.711 



0.388 





*3 


observed 


7t(D 


^predicted 


77(2) 


^predicted 


100 

200 

400 

800 

1600 

3200 

6400 

12800 

25600 


19.623 


JL 


13.049 

3.519 

1.058 

0.403 

0.219 

0.163 

0.144 

0.137 

0.134 


0.191 


JL 


B* - 97*A 

A 

"too ■ ^5° 

(L/R) 2 /Pe A - 436.2 

(1) Inter- Intraparticle Adsorption Model 

(2) Xntexpartide Equilibrium Adsorption Model 
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TABLE XV 


MOMENT ANALYSIS i*ND PARAMETRIC 
STUDY - ACETONE 100°C. 


»R 0 




„ ( 2 ) 

observed 

n 

^predicted 

T., predicted 

0.029 

173.29 

158.69 

156.49 ' 




_ <U 

( 2 ) 

m hj 

^ 2 , observed 

^ 2 , predicted 

^ 2 , predicted 

100 

815.67 

977-55 

437.28 

200 



723.1*1 



400 



686 . 3 I* 



800 



517.80 



1600 



483.53 



3200 



1 * 66 . 1*0 



61*00 



457.83 



12800 

25600 


/ 

453.55 

451>1 


f 

N RU- 


t 

_ ( 1 ) 

(2) 

^3, observed . 

^ 3 >predicted 

^ 3 >predicted 

100 

25404.0 

23392.2 

19499.2 

200 



20454.8 



1*00 



19745.3 



800 



19555.4 


. 

1600 



19531.7 



3200 



19485.1 



61*00 



19480.4 



12800 


f. 

19477.2 


f 

25600 

19476.2 




Pe A » 86 O 9 

"too • 88960 

(l/R) 2 /P« A - 3 * 8-2 



Inter- Intraparticle Adsorption Model 
Interparticle Equilibrium Adsorption Model 
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also desirable that the solution obtained by use of the finite terminal 
boundary condition possess the same characteristics at the column out- 
let (z = 1) as the infinite column solution, using Equation 4 as the 
terminal boundary condition. 

A possible alternate boundary condition sometimes encountered in 
similar diffusions problem is the zero-derivative condition 

<) y(z Q , 9)/ <^z = 0; z Q ^ 1 ; 9 > 0 (5) 


To develop guidelines on the use of Equation 5, two cases were compared 
via moment analysis: 

Case I Case II 


Equation 1, R = 0 Equation 1, R - 0 

Equations 2, 3, 4 Equations 2, 3, 5 

Specifically, parametric studies were made by moment analysis to 
determine 

1) errors associated with letting zq be 1 as a function 
Pe and 

2) the smallest z^, as a function of Pe that could be 
used to approximate infinite column behavior at z of 1. 

These results are summarized in Tab '1. Thus criteria have been 
developed for accurate use of a finic terminal boundary condition 
which can be accommodated in a numerical procedure. 

Considerable effort has been spent in evaluating several numerical 
methods for solving the partial differential equation set of Figure 91. 
The implicit Crank-Nicolson method, which is based on finite difference 
representation of the derivatives, appears impractical for this pulse 
propagation problem because of restrictions in space and time discreti- 
zations to maintain numerical stability. This appears primarily caused 
by the large Peclet number values (1000 to 10,000) encountered in these 
systems. 

The method of weighted residuals (MfR), Ref. 44, offers promise. 

In this method, the exact solution is approximated by a trial solution 
of the form: 



ii <*> ♦ f o <*) 


yJ'yOO 
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TABLE XVI. CRITERIA FOR USING FINITE BOUNDARY CONDITION 


Case I: infinite Boundary condition 

Case II: finite "boundary condition 



Safe z Q 

Safe Zq 

Pe 

<“ 1,1 = “l.n 1 * 

fi*2,I = U 2fll 

2 

9.791 

11.768 

4 

5.254 

6.021 

8 

3.043 . 

3.328 

i6 

1.978 

2.073 

32 

1.467 

1.490 

64 

1.223 

1.222 

128 

1.106 

' 1.044 

254 

1.050 

1.002 

512 

1.024 

1.0005 

1024 

1.011 

1.0001 

2048 

1.005 

1.0000 

4096 

1.002 

1.0000 

8192 

1.001 

1.0000 


I 


* “1,1 - “i,n - 10 ‘ 8 

** “2,1 * “2,h - “■* 
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where £^(z) are a set of trial functions and (a) satisfies 

the boundary conditions. If the trial solution is substituted into 
the partial differential equation, a residual, R, is formed. In MWR 
the weighted integrals of the residual are set to zero: 


J w. Rd V » 0 , j « 1, 2, ... N 

V J 

where wj is a weighting function. In the orthogonal collocation 
version of the method, the weighting functions are families of ortho- 
gonal polynomials. This method has been shown applicable to similar 
diffusion-type problems. Ref. 44. The application of orthogonal 
collocation to the system models under consideration is the present 
area of endeavor. It is anticipated that a technical report dealing 
with this subtask from July, 1972 through August, 1973 will be issued 
during September, 1973. 
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V. Educational Considerations 

In addition to the technical goals outlined earlier, this project 
has the objective of promt ing the education of engineering students in di- 
rections and dimensions not normally encountered in formal programs of study. 
In brief, the project activity provides real problems of substance whose 
solutions are obtained not by the application of expedient simplifications 
to make them manageable but rather according to the technical needs that 
realistically apply to the situation. Accordingly, the students learn to 
accommodate to trade-offs between competing values and to work with boundary 
conditions and constraints originating with the factual situations and which 
may be compounded by the interfacing of tasks. In this type of environment, 
the student perceives that his role in professional practice will normally 
involve a significant amount of interaction *7ith other individuals and that 
his work cannot proceed independently and without consideration of and im- 
pingement by the work of others. Furthermore, he is forced to obtain, under- 
stand and utilize knowledge which is on occasion far removed from his own 
speciality field and to undertake research when necessary to obtain the re- 
quired information. Although the very nature of the problems and the faculty 
perspective emphasize the relevance and importance of the technical goals, 
the periodic visits of Dr. Garrett Paine of Jet Propulsion, and on occasion 
of other NASA representatives, to review progress reinforce these concepts 
dramatically. 

From an educational point of view, the project has proven to be an 
unqualified success. Large numbers of students (approximately 100 in five 
years) have had the unique experience cf a real-life involvement in engineer- 
ing problems within the context of their formal education. Although the 
majority of participating students have been candidates for the Master of 
Engineering degree, a significant number have pursued the doctoral degree 
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and a substantial number of undergraduates have been involved. While a 
fraction of the students receive financial support from the project, the 
majority have elected to participate in order to receive academic cicdit. 

The latter students are self-supporting, hold fellowships or are retained 
as teaching assistants. The fact that so many of them have elected to work 
on the project to gain academic credit indicates the interest which the 
project has generated and the educational value that the participating 
students perceive to be related with this activity. 

During the past year from July 1, 1972 to June 30, 1973, a total 
of 34 students have been involved in the project. Of this group, 3 are 
following a doctoral program, 18 are, or were, candidates for a master's 
degree and the remaining 13 were undergraduates. A number of the under- 
graduates intend to continue their involvement in the project. 

This substantial interest by the students indicates that the NASA 
support Is contributing significantly and meaningfully to their education. 
Furthermore, the effort to resolve crucial problems related to an unmanned 
mission to Mars has been amplified by the technical contributions of stu- 
dents who are receiving no or very littlu financial support from the project. 
Thus, this NASA project has the unique benefits of providing technical 
support to national goals in space exploration and of enhancing engineering 


education at Rensselaer 


